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The best laboratory constraints on strongly coupled chameleon helds come not from tests of gravity 
' per se but from precision measurements of the Casimir force. The chameleonic force between two 

I nearby bodies is more akin to a Casimir-like force than a gravitational one: The chameleon force 

, behaves as an inverse power of the distance of separation between the surfaces of two bodies, just 

' as the Casimir force does. Additionally, experimental tests of gravity often employ a thin metallic 

' sheet to shield electrostatic forces, however this sheet mask any detectable signal due to the presence 

of a strongly coupled chameleon field. As a result of this shielding, experiments that are designed to 



specifically test the behaviour of gravity are often unable to place any constraint on chameleon fields 
with a strong coupling to matter. Casimir force measurements do not employ a physical electrostatic 
CO , shield and as such are able to put tighter constraints on the properties of chameleons fields with 

a strong matter coupling than tests of gravity. Motivated by this, we perform a full investigation 
on the possibility of testing chameleon model with both present and future Casimir experiments. 
^ ^ We find that present days measurements are not able to detect the chameleon. Ifowever, future 

^ ^1 experiments have a strong possibility of detecting or rule out a whole class of chameleon models. 
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One of the most common predictions made by modern theories for physics beyond the standard model is the 
existence of light scalar fields. It is usually the case that these fields couple to matter and hence mediate a new (or 



> 
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\ 'fifth') force between bodies. To date, however, no such new force has been detected, despite numerous experimental 
■ attempts to do so [1] . Any force associated with light scalar fields must therefore be considerably weaker than gravity 
over these scales, and under the conditions, that have so far been probed experimentally. This imposes a strong 
constraint on the properties of any new scalar fields; they must either interact with matter much more weakly than 
gravity does, or they must be sufficiently massive in the laboratory so as to have remained undetected. If the mass, 
, 7710, of the scalar field is a constant then one must require that hc/m^ ^0.1 mm if the field is to couple to matter with 
a strength equal to that of gravity. The bounds on fields whose interactions with matter have a super-gravitational 
strength are even tighter [1] . 

It has recently be shown, however, that the most stringent experimental limits on the properties of light scalar 
5-H ' fields can be exponentially relaxed if the scalar field theory in question possesses a chameleon mechanism [2, 3]. The 
chameleon mechanism provides a way to suppress the forces mediated by the scalar fields via non-linear field self- 
interactions. A direct result of these self-interactions is that the mass of the field is no longer fixed but depends on, 
amongst other things, the ambient density of matter. The properties of these scalar fields therefore change depending 
on the environment; it is for this reason that such fields have been dubbed chameleon fields. Importantly, Chameleon 
fields could potentially also be responsible for the observed late-time acceleration of the Universe [4, 5]. If this does 
indeed turn out to be the case, it raises the exciting prospect of being able to directly detect, probe and potentially 
even manipulate dark energy under controlled laboratory conditions. The properties of chameleon field theories are 
constrained by experimental tests of gravity, however, as a result of their chameleonic behaviour, theories in which the 
fields and matter interact with at least gravitational strength are not currently ruled out [2, 3]. Indeed, laboratory- 
based gravitational tests alone cannot even place an upper bound strength of chameleonic interactions with matter [3] . 
It was recently shown that some strongly-coupled (i.e. compared to gravity) chameleon theories predict alterations 
to the way in which light propagates through the vacuum in the presence of a magnetic field [6, 7]; the resultant 
birefringence and dichroism could be detected by laboratory searches for axion-like-particles e.g PVLAS, Q&A and 
BMV [9]. 

In Ref. [3] it was shown that the best laboratory constraints on strongly coupled chameleon fields come not from 
tests of gravity per se but from precision measurements of the Casimir force. In some ways this is not surprising. As 
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we shall sec, the chameleonic force between two nearby bodies is, in many ways, more akin to a Casimir-likc force 
than a gravitational one. Much like the Casimir force, the chameleonic force generally depends only very weakly on 
the composition and density of the test masses and in one class of theories, the chameleon force behaves as an inverse 
power of the distance of separation between the surfaces of two bodies. Additionally, unlike gravitational forces, the 
chameleonic force can be shielded [3] . 

Experimental tests of gravity often employ a thin metallic sheet to shield electrostatic forces, however this sheet 
was also shown in Ref. [3] to mask any detectable signal due to the presence of strongly coupled chameleon fields. 
This is because, in such theories, the shield develops what is known as a thin-shell. This means that the range of the 
chameleon field, X^j, = hc/m^, inside the metallic sheet is much smaller than the thickness of the sheet, rfghioid- In 
experimental tests of gravity, one measures the force or torque on one test mass ('the detector') due to the movement 
or rotation of another ('the attractor'). The electrostatic shield sits between the two. The shield is held fixed relative 
to the the detector and is uniform. As a result, residual forces due to the shield itself do not result in any detectable 
effect. In chameleon theories, the electrostatic shield attenuates the chameleonic force (or torque) due to the attractor 
by a factor of exp(— m^dshicid)- If m^dshieid ^ 1, the electrostatic shield therefore acts as a near perfect shield of 
the chameleonic force due to the attractor. Since is larger for strongly coupled fields than it is in more weakly 
interacting ones, experiments that are designed to specifically test the behaviour of gravity are often unable to place 
any constraint on chameleon fields with a strong coupling to matter. Casimir force measurements, on the other hand, 
do not employ a physical electrostatic shield and as such arc able to put tighter constraints on the properties of 
chameleons fields with a strong matter coupling than tests of gravity. 

A preliminary analysis of the constraints on chameleon fields provided by Casimir force measurements was made 
in Ref. [3]. In this paper we refine, extend and generalize this earlier study. Our primary aim is to extract the 
bounds that measurements of the Casimir force currently place on chameleon theories and to make predictions for 
what near future Casimir experiments will be able to detect. We shall see that there is a very real prospect that 
the next generation of Casimir force experiments will be able to detect or rule out most chameleonic models of dark 
energy. 

This paper is organized as follows: In Section II we introduce the chameleon model in greater detail as well as the 
concept of a thin-shell. When dealing with gravitational tests, it is generally the case that if the test masses have thin- 
shells, then all detectable effects due to chameleon fields are exponentially attenuated. In Casimir force experiments, 
however, the opposite is true. This is because when very small separations are used, the gradient of chameleon force 
is largest for thin-shelled test bodies, and Casimir tests are generally most sensitive not to the magnitude of any 
new forces but to their gradients. In Section III we present the conditions that must be satisfied for a test body 
to have a thin-shell. In Section IV we derive the form of the chameleonic force between two nearby bodies such 
as those used to measure the Casimir force. These results are applied in Section V to predict the extra force that 
should be detected by Casimir force measurements if chameleon fields exist. We also consider to what extent current 
experiments constrain two of the simplest and most widely studied classes of chameleon theories. In the penultimate 
section, we then consider the extent to which planned future experiments will be able to extend the constraints on 
chameleon theories, and identify two proposed tests that have the sensitivity to detect or rule out most chameleon 
theories in which the chameleon potential is associated with dark energy. We conclude in VII with a discussion of our 
results. 



As was mentioned above, chameleon theories are essentially scalar field theories with a self-interaction potential 
and a coupling to matter; they are specified by the action 



where (p is the chameleon field, Sm is the matter action and Vm are the matter fields; V{(f)) is the self-interaction 
potential. 

The strength of the interaction between and the matter fields is determined by the one or more mass scales Mj. 
In general, we expect different particle species to couple with different strengths to the chameleon field i.e. a different 

Mi for each V^m- Such a differential coupling generally leads to violations of the weak equivalence principle (WEP 
hereafter). Constraints on any WEP violation are very tight [1]. Importantly though, it has been shown that V^((/') 
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can be chosen so that any violations of WEP are too sniaU to be have been detected thus far [2, 3]. Even though the 
Mi are gcneraUy different for different species, if ^ 0, we expect Mj ~ 0{M) where M being some mass scale 
associated with the theory. In this paper we are concerned with those signatures of chameleon theories that could be 
detected through measurements of the Casimir force. Since these measurements place bounds on the magnitude (or 
gradient) of close range forces rather than on any violation of WEP, and since also all that matters in this context is 
the coupling of the chameleon field to atoms rather than any more exotic form of matter, allowing for different Mi is 
an usually an unnecessary complication. Henceforth, we assume a universal coupling M^ = M for all i and take the 
matter fields to be non-relativistic. The scalar field, then obeys: 

Ucj> = V\^) + ^, (3) 

where p is the background density of matter. The coupling to matter implies that particle masses in the Einstein 
frame depend on the value of (j) 

m{(j)) = e'f'/^mo (4) 
where mo = const is the bare mass. We parametrize the strength of the chameleon to matter coupling by (3 where 



and Mpi = l/VSwG w 2.4 x 10^^ GeV. On microscopic scales (and over sufficiently short distances), the chameleon 
force between two particles is then 2/3^ times the strength of their mutual gravitational attraction. 

If the mass, = ■\/V'"(^), of is a constant then one must either require that > 1 meV or /3 <C 1 for such 
a theory not to have been already ruled out by experimental tests of gravity [1]. If, however, the mass of the scalar 
field grows with the background density of matter, then a much wider range of scenarios is possible [2-4]. In high 
density regions can then be large enough so as to satisfy the constraints coming from tests of gravity. At the 
same time, the mass of the field can be small enough in low density regions to produce detectable and potentially 
important alterations to standard physical laws. Scalar fields that have this property are said to be Chameleon fields. 
Assuming dlnm(0)/d0 > as it is above, a scalar field theory possesses a chameleon mechanism if, for some range 
of (j), the self-interaction potential, V{^), has the following properties: 

V'{(l)) < 0, V" > 0, V"'{^) < 0, (6) 

where V = dV/d(p. Whether or not the chameleon mechanism is both active and strong enough to evade current 
experimental constraints depends partially on the details of the theory, i.e. V{(j)) and M, and partially on the initial 
conditions (see Refs. [2-4] for a more detailed discussion). For exponential matter couplings and a potential of the 
form: 

y(0) = A^ exp(A"/0") ^At(^l + (7) 

the chameleon mechanism can in principle hide the field such that there is no conflict with current laboratory, solar 
system or cosmological experiments and observations [2, 4]. Importantly, for a large range of values of A, the chameleon 
mechanism is strong enough in such theories to allow even strongly coupled theories with M <C M pi to have remained 
undetected [3]. The first term in V{(l)) corresponds to an effective cosmological constant whilst the second term is a 
Ratra-Peebles inverse power law potential. If one assumes that ^ is additionally responsible for late-time acceleration 
of the universe then one must require Ao ~ (2.4 ±0.1) x 10""'^^ GeV. In the simplest theories, A ^ 0{Aq) so that 
there is only energy scale in the potential; it is arguable that this represents the most natural scenario. The smallness 
of Aq, means that, as a dark energy candidate. Chameleon theory do not solve either the naturalless problem or 
the coincidence problem. However, although it would certainly be desirable to have a model which solved both of 
these problems, one cannot exclude the possibility that the acceleration of the Universe is the first sign of some new 
physics associated with an 0(Ao) energy scale. If this is truly the case then one must look for new ways in which to 
probe physics at this low energy scale. As we show in this paper, the use of Casimir force experiments to search for 
chameleon fields is once such probe. 

Throughout the rest of this paper, it is our aim to remain as general as possible and assume as little about the 
precise form of V{(j)) as is necessary. However, when we come to more detailed discussions and make specific numerical 
predictions, it will be necessary to chose a particular form for V{(j}). In these situations we assume that V{<f)) has 
either has the following form: 

/ A" 

y(<^) = AMi + — 
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or 

y(</.)=A^exp(A"/.^"). 

We do this not because these forms of V are in any way preferred or to be expected, but merely as they have been the 
most widely studied in the literature and, in the case of the power-law potential, because is the simplest with which 
to perform analytical calculations. The power-law form is also useful as an example as it displays, for different values 
of the n, many the features that we expect to see in more general chameleon theories. 

The evolution of the chameleon field in the presence of ambient matter with density pmatter is determined by the 
effective potential: 

V,tt{4>) = V{4>)+p^,ucre'^/^ (8) 

Limits on any variation of the fundamental constants of Nature mean that, under the conditions that arc accessible 
in the laboratory, we must have (p/M <C 1 [3]. Henceforth we therefore take exp(0/M) « 1 + (j)/M. Even in theories 
where V{(j)) has no minimum of its own (e.g. where it has a runaway form), the conditions given by Eq. (6) on 
ensure that the effective potential has minimum a,t (f) = 0min(Pmattcr) where 

Kff('/'mi„) = = V'{4>min) + (9) 



B. Thin-shells 

In chameleon field theories, macroscopic bodies may develop what has been called a 'thin-shell'. Generally speaking, 
a body of density pc is said to have a thin-shell if, deep inside that body, cj) is at, or lies very close to, the minimum 
of its effective potential (where (j) = 4>c ^ (pminipc) say). We take the density of matter outside the body to be pb] far 
outside the body (f> (pi, = (pminipb)- 

Thin-shelled bodies are so-called because, for such bodies, almost all of the change in (p (from </>;, to 0c) occurs in 
a thin region near the surface of the body. The thickness of the part of this thin region that lies inside the body is 
generally » 0{l/mc} where m,; = m^{4>c)- If the body has thickness R then, a necessary (but not sufficient) condition 
for a body to have a thin-shell is nicR » 1. 

As a rule of thumb, larger bodies tend to have thin-shells whereas smaller bodies do not. Precisely what is meant 
by 'larger' and 'smaller', however, depends on the details of the theory. We discuss this further in Section III below. If 
mft = 7710 (0t) is the mass of the chameleon in the background then the chameleon force between two non-thin-shelled 
bodies, separated by a distance r, is 2/3^6"™''' times as strong as their mutual gravitational attraction. Importantly, 
the force between two thin-shelled bodies is much weaker [2]. Moreover, it has been shown that the chameleonic force 
between two thin-shelled bodies is, to leading order, independent of the strength, /?, with which the chameleon field 
couples to either body [3]. In a body with a thin-shell, it is as if the chameleon field and the resultant force only 
interact with and act on the matter that is in the thin-shell region near the body's surface. 

III. THIN-SHELL IN CASIMIR FORCE EXPERIMENTS 

Before we can consider the form or magnitude of chameleonic force wc need to know whether or not the test 
masses used to measure the Casimir force are predicted to have thin-shells. In subsection III A we state the thin-shell 
conditions for an isolated spherical body in the context of general chameleon theory, which are themselves derived 
in Appendix A. Wc then consider whether these conditions hold for the test masses used in those Casimir force 
measurements that have been conducted thus far in subsection IIIB. 



A. Thin-Shell Conditions 

In general, 'larger' bodies have thin-shells whereas 'smaller' ones do not. How small is 'small', however, generally 
depends on the details of the theory. The test-masses used in Casimir force experiments come in a number of different 

shapes and sizes. Some experiments used relatively small test masses, with typical length scales of 0(10^ whilst 
others, perhaps most notably that performed by Lamoreaux in 1997 [10], used relatively large test masses with length 
scales of 1 — 10 cm. 

The condition that must be satisfied for an isolated spherical body to have a thin-shell was first derived, for V oc (^~" 
with n > potentials, in Ref. [2]. In Ref. [3], the thin-shell conditions for such potentials were re-derived (via a 
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different method) and extended to theories with n < —4. The thin-shell conditions for theories with n < —4, n = —4 
and n > were found to be qualitatively different. 

In Appendix A, we derive the thin-shell condition for general V{(j)). We consider an isolated spherical body with 
density Pc, radius i? in a background with density pb- We define (pb by V'{(j)b) = —pb/M and 0c by V(?^c) = —Pc/M. 
We also define rub = m^{(j)b) and nic = m^{(j)c)- If nibR ^ 1 then the body will always have a thin-shell of some 
description as almost all variation in cj) will take place in a thin region (of thickness at most ~ 1/ mb) near the surface 
of the body. However, if rric « mb then this thin-shell would be linear, i.e. we would see almost the same behaviour 
if we considered a Yukawa theory with mass nib (for which the field equations would be linear). In the cases we 
consider, however, pc 3> Pb and so necessarily nic ^ m^. If rricR is large enough then, whatever the value of rribR, a 
body may have a non-linear thin-shell. Non-linear thin-shells are associated with the dominance, near the surface of 
the body, of non-linear terms in the field equations. Such behaviour would not occur in theories where (f) has only a 
Yukawa coupling to matter. This non-linear behaviour is key in allowing chameleon theories to evade the stringent 
experimental constraints on the coupling of a scalar field to matter that exist for Yukawa theories. In Appendix A 
we find that a necessary and sufficient condition for a non-linear thin-shell, in a general chameleon theory, is: 



C = 



(pc - Pb)f{mbR)R^ 



> 



2M[mlR^ + rubR + 1] 
where this defines C. An equivalent statement of this condition is: 



{Pc- Pb)i'i-- fjmbR)) 



mlR^ 



> 



2Mml {(fb - <t>o) 2(1 - f{mbR)) 



mlR^ + rribR + 1 ~ (pc - Pb)f{mbR) 



f{mbR) 



> 2. 



(10) 



(11) 



We have defined: 



f{mbR) = 26-"'"^ coshinibR) 



rribR 



mlR^ 



tanh mfei? 
rubR 



As lUbR 0, f{mbR) 2/3 and as mbR. oo, f{mbR) 
Far from a spherical body with a non-linear thin-shell ( 



1. 

has the form: 



<Pb - 



-'thin 



where 



CtMn{l + mbR + mlR^) 
V'{cl>b)-V'{cj>b-Cti,in) 



R^ 
2 ■ 



We note that Cthin = Cthin(-R, 4>b) and is therefore independent of Pc/M. It follows that, in all chameleon theories, far 
from a body with a non-linear thin-shell (/) is independent of Pc/M i.e. it is independent of the strength with which 
the chameleon field couples to the matter in the body. This, more than anything else, is what makes it so difficult for 
experimental tests of gravity to place a lower bound on M. 

In many cases, it is only necessary to consider the following sufficient condition for a thin-shell: 

C>,d)b- (be. 



where C is given by Eq. (10). 



B. Applying the Thin-Shell Conditions 

We now consider whether or not the test masses used in experimental measurements of the Casimir force are 
generally predicted to have non-linear thin-shells. The above thin-shell conditions are valid for isolated, spherical 
bodies. Generally, however, at least one of the test masses used in Casimir force measurements is neither isolated nor 
spherical. Isolated in this case means that there is enough space between the body in consideration and any other 
bodies for, in all directions, m,f,{4') to have decreased to be about rub before any other body is encountered. 

Generally speaking, — > nib over a distance of about l/mb- We therefore take an isolated body to be one outside 
which there is a region of thickness at least > l/rub in which pi^ pb- 

For isolated non-spherical bodies, such as rectangular plates with volume V and longest dimension 2D, the above 
thin-shell conditions still apply (to a good approximation) provided one replaces R by ^JW/AttD (i.e. one should 
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replace by the volume divided by the longest distance from the centre of mass of the body and its surface). If 
two plates with volumes Vi and V2 and longest dimensions 2Di and 2D2 respectively are placed a distance d apart, 
with d > l/rrib, Eq. (10) gives the thin-condition for each plate with R replaced by ^3Vi/4nDi for i = 1,2. If, 
however, d <C 1/mb, then plates are not isolated and they effectively count as one mass for the purposes of applying 
the thin-shell conditions. Provided l/rric is small compared to the smallest dimension of each plate, the thin-shell 
condition for both plates is then given by Eq. (10) but with R = y^314ot/47i^tot where 14ot is the total volume of the 
plates (excluding the space in between them) and Dtot is the half longest dimension of the two plates when considered 
as a single object. 

The overall geometry of the set-ups used for Casimir force measurement is generally quite complicated . Even in 
experiments where the test masses are themselves relatively small and thin, the apparatus that surrounds them is not. 
Furthermore, the test masses are generally not isolated in the sense defined above. This complicates the application 
of the thin-shell conditions, and generally it can only be done thoroughly within the context of a specific chameleon 
theory or class. 

For definiteness and as an example we consider theories where V{(j)) has a Ratra-Pebbles form V{(j)) = Ag(l -|- 
{A/([))"), where n > 0. We take Aq = 2.4 x 10~^eV so that the constant term in the potential is responsible for the 
late time acceleration of the Universe and specifically consider theories where A w Aq. In Casimir force experiments 
the test masses are much denser than the laboratory vacuum in which they sit and so wc take ^ pi,; this implies 
that (j)c -C (j)b since n > 0. Given these considerations, the thin-shell condition, Eq. (10), simplifies to: 

^ _ PcR^finibR) ^ 



2M{l + mbR + mlR'^) 

Since pc ^ Pb if Wfc-R > 1 then this condition is automatically satisfied. We therefore restrict our attention to those 
cases where rribR <C 1. For pc » Pb and rribR <C 1, the thin-shell condition is: 



GMbody 



> 



R ^ 2/32M' 

where Mbody is the mass of the body and /? = Mpi/M. Note that such a simplification of the thin-shell condition will 
generally occur for all theories where {(f) can be shifted so that) ^{p) — > as p = —MV'{(f){p)) — > 00. 

We take the pressure of the laboratory vacuum to be p x 10~^torr, p ^ 0{1) or greater for all Casimir force 
measurements made to date. We then find that: 

= f3 r^+lp n + 



2/32M 
where 

Bn = 4.9 X 10"^^ (5.1n x 10^°) ^ . 

The largest value of Bn occurs for n 0.048 at which Bn ~ 4.5 x 10~^^ and Bn decreases very quickly to 4.9 x 10"^"'^ 
as n ^ 00. Additionally, as is discussed more fully in Refs. [2, 3], one must be aware that the smallest value that nib 
can take in laboratory vacuum which has a smallest length scale Lvac is 0(1/Lvac)- If a typical value of ivac = 1 m, 
we therefore have 



^bf^P M < jj-^ 4.9 X 10-31 (l.5n(n-|- 1) x 10^) 



< 1.5 X 10"^^. 



For bodies with p^ > Pgiass ~ 3 gem ^, we have: 



6.7x10-27^ ^ 



R \R cm, 

where V is the volume of the body. If we take V to be the volume of the smallest isolated system associated with 
either of the test masses, and R = I?iong the longest distance from the surface of this system to its centre of mass, 
then the test bodies will certainly have thin-shells if: 



/ V 1 
-Reff = \ I — 7^ cm. 

V i'long ~ 2^ 

For this choice of potential and \n\ ^ C(l), the Eot-Wash experiment [14] requires that if /3 > 10"^ and A = 
2.3 X 10-3 eV then (3 must be larger than 10^ [3]. For V = Ag(l -|- A" /</>") with A = Aq = 2.4 x IQ-^eV, ||n|| ~ 0(1), 
we have checked that the thin-shell condition certainly holds for the test masses used in the Casimir force experiments 
reported in Refs. [10, 15-19] provided /3 > 10"^. In other words, it holds for most strongly coupled chameleon theories 
that are not already ruled out by tests of gravity such as the Eot-Wash experiment. 
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C. Discussion 

We found above for V = A^(l + A" /(f)), \n\ ~ 0(1) and A = Aq = 2.4 x 10^ eV, the test masses used in all 

Casimir force experiments conducted to date are predicted to have thin-shells in all theories with /? ^ 10'^. In some 
experiments, where particularly large test masses are used, the test masses are also predicted to have thin-shells for 
0(1) values of /3. 

1{V = A4/((A/(/))"), n > 0, A « Ao = 2.4 X lO'^ eV where /" > and / is normalized so that /' = 1 (e.g. 
V = Ag exp(A"/(/)"), then the potential is always steeper than the Ratra-Peebles form considered above and as such 
the thin-shell conditions are less stringent. 

In the next section we calculate the chameleonic force between two nearby bodies under the assumption that they 
have thin-shells. In the absence of thin-shells, the chameleon field behaves in the same way as a Yukawa field, and the 
constraints on any Yukawa coupling to matter derived from Casimir force measurements can be directly applied to 
chameleon theories. There is, therefore, nothing new to say about the non thin-shelled case, and we do not consider 
it further. 

We find below that the gradient in the chameleonic force between two nearby thin-shelled bodies is gc;nerally much 
steeper than it would be if there were no thin-shells present. Casimir force experiments generally measure gradients 
in forces (changes in forces between two separations). So they are generally more sensitive to relatively small quickly 
varying forces with a steep gradient than they are to large but nearly constant ones. The presence of thin-shelled test 
masses is therefore an aide rather than a hindrance to the detection of chameleon fields via Casimir force measurements. 
The opposite is generally true of gravitational tests [2, 3]. The stronger the matter coupling, the more likely it is that 
a given body has a thin-shell. Experiments designed along the lines of Casimir force measurements are therefore far 
better suited to the search for strongly coupled chameleon fields than tests that are specifically designed to search for 
forces with gravitational (or sub-gravitational) strength. 

IV. THE CHAMELEONIC 'CASIMIR' FORCE 

In this section we calculate the Casimir-like force between two nearby bodies due to their interaction with a 
chameleon field. The form of both the Casimir force and the chameleonic force are highly dependent on the geometry 
of the experiment [3, 11, 12]. The form of these forces is most easily calculated when the geometry is that of two 
parallel plates. Making accurate measurements of forces using this set-up is, however, notoriously diflacult as it 
requires that the plates be both very smooth and held parallel to a high precision. For this reason, most experiments 
conducted to date have measured the Casimir force between a plate and a sphere rather than between two plates. 
Presently the highest precision measurements have been made using the sphere-plate geometry. By measuring the 
gradient of such a force between a sphere and a plate it is possible to determine the force between two parallel plates. 
In Section IV A we calculate the chameleon force between two parallel plates, and in Section IV B we calculate the 
chameleon force for the sphere-plate geometry. 

A. Parallel Plates Geometry 

The parallel plate geometry is the easiest to study analytically. For simplicity we take both plates to have the same 

composition. We shall see that, provided the plates both have thin shells and are much denser than their environment, 
the chameleon force is in largely independent of composition of either plate. We take the plates to have thin-shells 
and to be separated by a distance d. 

We define a; = to be the point midway between the two plates; the surfaces of the plates are then at a; = ±d/2. 
In —d/2 < X < d/2 we have: 

g = v'{4>) - v'iM, 

where we have defined V'{(j)b) = —pb/M with pj, the ambient density of matter outside the plates. Inside the plates 
we have: 

where V'{(f)c) = —pc/M and pc is the density of the plates. Deep inside either plate ^ const w (pc and dcp/dx = 
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at a; = by symmetry. We use the shorthand (?io = (p{x = 0). Integrating both of these equations once gives: 

(g)' = 2(V(0)-y(0o)-V'((/.b)((/)-0o)) -d/2<x<d/2, (12) 

(ai)' = 2 iv{<l^) - V(</>c) - v'{<p,){^ - > d/4. (13) 

We define (f)s = 4'{x = ±d/2), so that is the value of (j) on the surface of the plates. By matching the above 
equations are x = ±d/2, we arrive at: 

, _ V{ct>c) - r (0c)0c - V{M + V'{<t>b)^o ,,,, 
y'(0.) - y'(^e) • ^'^^ 

If one of the plates were to be removed then (pa on the surface of the remaining plate (= (pao say) would be given by 
Eq. (14) but with 0o ^ 4>b- The perturbation, 5<j)s = 4>s — 4>sO, in (ps due to presence of the second plate is therefore: 

^ v{(Pb)-v{M-v'i(l>b){<l>b-^o) 

V'i(l)b)-V'{(bc) 

Deep inside either plate the perturbation, d(j), in (j) due to the presence of the second plate is exponentially attenuated. 
This is because the chameleon mass inside either plate nic = m^{(pc) is, by the thin-shell conditions, large compared 
to the thickness of the plate. 

The attractive force per unit area, F^/A, on one plate due to the other is given by: 

dx— — « y ((/.e)(5?!>«, (15) 



A 7^/2 M dx 

where D is the plate thickness and we have used Pc/M = —V\(j)c)- Taking pc ^ pb we then find: 

^ = v{M - v{M + v'ihm - M < v{M - v{ci>b). (16) 

To leading order in Pb/Pc, F<i>/^ depends only on </>□ and (/>b. 

To a first approximation, we calculate (^o by linearizing the equation for (j) in —d/2 < x < d/2 about In 
-d/2<x< d/2 we find: 



^-^o=^^^sinh^(^). 



(17) 



"0 



Vq = V{(pa} and = V\(pb}- If, as is the case in theories with V = Af^f{{A/(f))^), n > and /' > 0, we 
expect that <p on the surface of the plates is very small compared to (po, then toq is given approximately by: 

sinh^f!!^)^-!!^. (18) 



V 4 ; 2{V^-V^) 
liV = A^(l + A"/?^") (and n > 0) when m^d » 1 this gives: 



mod « 4 sinh ^ 



(n + 1) 



This is a good approximation for n ^ 0{1), but it breaks down for larger values of n. More generally, we must 
calculate mp d by a more complica ted method that takes proper account of the non-linear nature of the potential. We 
define y = \/V'^^Vo^^\^{^^-'^ and 1/W{y) = (V^ — ¥'{(/))) > 0. Prom the chameleon field equation we then have: 



V2j^'^W{y)dy=^, (19) 



where t/g = y{(p = cpg). The above integral can then either be calculated numerically for a given V{4i) or, as is often 

more helpful, via an analytical approximation. We show below that for mribd ^ 1, mod ^ When ^ TOq, 

which therefore corresponds to rricd ^ 1, we have moysWa ^ 1. It can be checked that, W{y) always decreases faster 
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than l/y as y increases for y l/moWo- We can therefore approximate Eq. (19) by replacing j/g with oo as the 
upper hmit of the integral: 



W{y)dy = -, 



Since mo < rric, if rricd < 1 then we must have mo « TOc- In the mcd ^ 1 case, we proceed by defining 



where a subscript indicates that the quantity is evaluated for = ^g, and a subscript b means that it is evaluated 
for (j) = When mo '> rrib, which we shall see corresponds to m^d <C 1, k"^ is, for many choices of potentials, almost 
independent of 0o- The value of fc^ dictates the dynamics of the theory, and determines how (j)o depends on d. We 
evaluate Eq. (19) approximately in Appendix B. We find that if, as is often the case, 1/3 < /c^ < 2, then we can 
define 

rieff = (2 - _ 1) (21) 

leading to 



rieff V 2 ' 2 rieflF 

where B{-, •) is the Beta function. 1/3 < fc^ < 2 implies that Ueff > or Ueff < —5/2. This approximation becomes 
exact when V — A"^ + A'*(A/c/))" and mo 3> mt,. For such a potential Ueff — n, and the requirement that mo S> mi, 
implies m^d <C 2 (or all n). For very steep potentials, e.g. V = A'^ exp((A/0)") when (/)o <C A, fc^ « 1 when mo » m(,. 
This corresponds to n^jj — > oo. It is clear that mod ~ 0{1). 
If fc^ > 2 then we find in Appendix B that: 



mod : 



^3/2 



2^/2(fc2 -2)(l/2) 

where J_i/4(-) and y_i/4(-) are Bessel functions. For small 4(A;^ — 2) this gives: 

mod « (1 - 3(/e2 - 2)/8) , 

and if fc^ 2 we have: 



(23) 



mod i 



g(i.i) 5.24 
V2fc ~ Vfc ' 



The 1/3 < /c^ < 2 and fc^ > 2 approximations for mod are continuous at fc^ = 2. 

In Appendix B, we show that when fc^ < 1/3 the analytical approximation used to evaluate 0o(d) for fc^ > 1/3 
breaks down. When is small it is either because at least one of fcg = V^f^VQ/mf^ or 1 — VI/Vq is small. Theories 
with small k^ only exhibit very weak non-linear behaviour near 0o and so we do not consider them further. 

It is important to know how F(i,/A behaves as fc^ — > because as d — > oo we have (I)q ^ (j)b. In Appendix B we find 
that for mftd 5 we have: 

mbd, 



and so 



ln(12) - hiiP), (24) 



72mge-'»'-^ 

-J ^^^^ 



We have now derived the expressions for Fff,{d)/A for different classes of theory and for different ranges of d. These 
expression generally consist of an exact expression for /A as a non-linear function of (po and an approximate implicit 

equation for (f)Q as a non-linear function of d. In these cases an explicit equation for F^/A as a function of d can 
only be found once V{4>) is specified. In the limit d ^ oo, it was possible to find an explicit expression for F^/A 
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Behaviour of Chameleonic Pressure for V = aJ(1+a"/(^"); n = 1 Behaviour of Chameleonic Pressure for V = aJ(1+a"/(^"); n = 4 




Separation of plates: m d Separation of plates: m d 



Behaviour of Chameleonic Pressure for V =A^(1+a"/4)"); n = -8 




10"^ 10° 10^ 10* to'' io' 



Separation of plates: m d 



FIG. 1: The dependence of the chameleonic pressure, F^/A, between two parallel plates on separation, d. We have taken 
V{ij}) — Ao(l + A"/(^") and fixed nic/mb = 10^. The three plots show the behaviour of F^/A for a theory with n = 1, n = A and 
n = —8. Each of these are respectively representative of theories with 0<n<2, n>2 and n < —4. Three types of behaviour 
are clearly visible in these plots. For d < , F^/A « Vc — Vt which is independent of d: this is the 'constant force behaviour'. 
For < d < , F^jA oc for some p. Theories with < n < 2 have 0<p<l. Ifn>2 then 1 < p < 2 and if 

n < —4 we have 2 < p < —4. This is the 'power-law behaviour'. Finally when d -C in^^ ^ F^/A oc exp(— mi,d), i.e. we have 
'exponential behaviour'. Note that in a standard Yukawa scalar field theory (where = const) one would have F^/A ~ const 
for d <C mT^ and an exponential drop-off for d > mT^; however there would be no region of power-law behaviour. 



as a function of d. We note that when d ^ iti~^, the leading chameleonic force is independent of rric, and hence 
also of the strength with which the chameleon field couples to the body. In general, the approximate expressions 
for 4>o{d) derived above should be seen as providing a good order of magnitude estimate for F^/A rather than an 
accurate numerical prediction. This said, in some cases the expressions found above are actually exact. Specifically, 
if V{(l)) = Aq(1 -I- A"/(/)") then Eq. (22) is exact (and rie// = n) in the hmit m^^ < d < m~^. Additionally, if 
V{(j)) — Aq exp(g((/)/A)) and (I)q is such that g"{4)oA)/g''^{(j)oA) <C 1 then Eq. (22) with P = 1 provides an excellent 
approximation when d When one only wishes to consider a specific form of V{(j)), numerically 

accurate predictions for F^/A can be made by performing the integral in Eq. (19) numerically. 

We now consider two specific examples. Firstly, if V{(j)) = Aq(1 + A"/(/)") and m^^ <C d <C m^^ then: 

^ = A4x„(A,rf)-^, (26) 
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where = Ag/A and 





(27) 



In all chameleon theories with power-law potentials, F^/A drops off as l/dP for some p in this regime. In theories 
with 0<n<2, 0<p<l;ifn>21<p<2 and if n < -4 we have 2 < p < 4. 
If d < then mo ~ rric and so: 



A 



where the last approximation holds if pc ^ Pb- In this regime, the chameleonic force is independent of d at leading 
order. Finally if d m^^ we have from Eq. (25): 

72n{n + l)Vbe-""''^ 
~A (n + 2)2 ■ 

In FIG. 1 we show the behaviour of the chameleonic pressure, F^/A, in all three regimes {d < m~^, <C d <C m^^ 
and d » ?ti^^) for chameleon theories with n = —8, n = 1 and n = 4. In all cases we have fixed mc/mb = 10^. The 
n = —8 plot is representative of theories with n < —4 and the n = 4 is representative of theories with n > 2. The 
n — 1 graph show an example for a theory with < n < 2. The three types of behaviour: constant force for d < w""'^, 
power-law for d ''^b^ exponential drop-off for d <C are clearly visible in these plots. The main 

difference in the behaviour of the force for the different values of n is the slope of F^/A in the power-law drop-off 
region. 

If V{(j}) = A^ exp(A"/(^") and again m^^ «: d «: m-\ then: 



n^l{n + 1) -f 2(,^o/A)" + n(0o/A)' 



(<^o/A) 



In 



For small (A/^)", we have rieflf = n and hence F^jA is given by Eq. (26); this limit corresponds to mo <C Ad so 
d ^ Ad~^. In the opposite limit when mo » A^ i.e. d <C A^^, we have instead nes ~ n^{A/(j))'^"/{n + 1) and so: 



mid^ 



2-k' 



1 -41n(2) 



"eff 



It follows that: 



n^h{Add)^d'^ 



nh{Add) 



where /i(Adrf) is a slowly varying function of A^d defined by: 



27r2 



n2(A2rf/A)2- 



(28) 



(29) 



The above expression for F^/A is valid provided that {n + 1) / nh{Add) -C 1. Note that, in all cases, F^{d)/A ^ 0{Aq) 
when w A/Aq. When d < m~^, we have in all cases that: 



A 



Vc-Vb + vMb - (/-c). 



When d <C Tfi^^ ^ Eq. (25) gives the behaviour of F^,! A. 

Note, for comparison, that the Casimir force per unit area between two parallel plates with separation d at zero 
temperature is: 
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B. Sphere-Plate Geometry- 



Calculating the chameleonic and Casimir forces is simplest for the parallel plate geometry. However, the most 
accurate measurements of the Casimir force have been made using a sphere and a plate. In this geometry, the Casimir 
force depends only on the radius of curvature, R, of the curved body, and the distance, d <C i?, between the surfaces 
of the two bodies at the point of least separation. 

We now calculate the chameleonic force between a sphere with radius R and a circular plate with total surface area 
A. d is defined to be smallest separation between these bodies. The z direction is defined to be perpendicular to the 
plate, and we take r to be a radial coordinate which measures the distance from the point of least separation in the 
plane of the plate. Both the plate and the sphere are assumed to have thin-shells. 

In the region between the two bodies the chameleon field satisfies: 



d</.2 d2<^ 1 d</. 

At r, the separation between the sphere and the plate in the 2:-direction is 

s{r)=d + Rll-(l-^ 



(31) 



We define (ppp{z,s) to be the value of the chameleon field in the parallel plate set-up for plates with separation s. 
When s <^ R, we can approximate (t>{z,r) in the sphere-plate geometry by (?iipp(z, s(r)). This approximation is valid 
so long as: 



r dr 



d2(^ 
dS' 



< 1. 



Now in the parallel plate set-up we found that: 

/■y(0(z);0o) 

^2 



where z = is the surface of the plate. We then have: 

-W{y) y 



d^ 
d^ 



W(0) W{Q)Jq y' du 



1 dW 



ds 



y 

72" 



Therefore if > 1: 



ds 



V2 



— -dy' 
dy' y' 



(32) 



Inserting the equation for s(r) we arrive at: 

dV 



d# = t 



d(/) _ v^i?2 ^{R + d- s)2 A(t> 



R + d- s{r) 



ds' 



It is clear then that the approximation (j) « (j>pp{r, s{r)) is good provided that d < s <^ R i.e. r, d <^ R. All 
sphere-plate Casimir measurements have d <C i?. 
Whenever (p{z, r) w <ppp{z, s{r)) we have: 



dA 



and for r <^ R, dA ~ 2'jTrdr « 2TrRds. The contribution to total force between the plate and the sphere from the 
points with r <C i? is: 

F^{r) « 2nR ^ ds' [V^(</'o(s(r'))) - V^(</'6) - V^'(.^6)(<^o(s(r')) - <l>b)] ■ 
Jo 
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We define rmax = min(-\/ A/tt, R) i.e. Tmax is the smaller of R and the radius of the circular plate. We also define 
Smax = s(rijiax); Sniax IS then the largest separation between the surfaces of the plate and the sphere. If dF^/dA drops 
off faster than 1/s for all s > s* for some s* ^ R, then the dominant contribution to the total force between the 
sphere and the plate comes from the region where s <^ R. To a very good approximation we therefore have: 

/■Smax 

« 27rR / ds' [V{Ms{r'))) - ^(<^6) - V'{<t>b){Msir')) - h)] ■ (33) 

Jd 

In some chameleon theories, however, dF^/dA drops off more slowly than 1/s for all s < 0{R) <C m^^. In these 
cases, Eq. (33) is no longer accurate. 

In all theories dF^°^{d)/dd w 2TTRdF(j,{d)/dA provided d R. By measuring the gradient of F^°*(d), it is therefore 
possible to extract the form of dF^{d)/dA. dF^{d)/dA is equal to the force per unit area between two parallel plates 
with distance of separation d. 

It is clear that the dependence of dF^/dA on s plays an important role. In particular, theories where it drops off 
more slower than 1/s behave differently from those where the drop off is faster. When toq 2> Wf, we have: 

^^V{<l>o)-V{ct>t,). 

We define Q{(j)o;V) = dln(y(0o) - F(<?!),,))/dln(l/s). If Q{Ms);V) > 1 for all s > s* (and < 1 otherwise), where 
s* <C min(smax, l/wfc), then the dominant contribution to the total force comes from points with separations « s*. In 
these cases Eq. (33) provides a good approximation to If no such s* exists but mbSmax ^ 1 then the dominant 
contribution to comes from points with separations w l/mfe; in these cases we may also use Eq. (33) to calculate 
If neither of these conditions hold, then dominant contribution to F^°^ comes from separations ^ ^^(sinax)- The 
assumption that <p ~ ^pp(^, s(r)) fails for s ~ C(si„ax) and it is particularly bad if Smax ^ 0{R). In these cases the 
chameleon field equations are too complicated to solve analytically. However, Eq. (33) is still expected to provide 
an order of magnitude estimate for F^^ . This is because the assumption that « 0pp(z, s(r)) only breaks down for 
s ~ C'(smax) but holds for all smaller values of s. We therefore, do not expect 0(z, Smax) to be very different from 
(^pp(i;, Smax) or to be very different from the form given by Eq. (33). 

liV = Ag(l + A"/(/)") then dF^is) / dA drops off more slowly than 1/s when m^s <C 1 for —2 < n < 2. Theories 
with —2 < n < are not valid chameleon theories. Theories where < n < 2 therefore make qualitatively different 
predictions for than do those where n > 2 or n < —4. 

In the sphere-plate geometry (with d R, m^"^ and nibd <C 1), the total chameleonic force for theories 

with n > 2 or n < —4 is given to a very good approximation by: 

Ff\d) « 27rA2Aii (^^) Kn {Addy^^ , (34) 

where Kn is given by Eq. (27) and, as above, A^ = Ag/A. 
In theories with < n < 2, however, we have: 

Fl°\d) « Fo(smax, mb) - 2TTKlhR ( (A^d)^ . (35) 

where i^o is independent of d and is calculated in Appendix C. If mfeSmax 1 then we are only able to find the order 
of magnitude of Fq: 

Fo~27rA^Ai?(^^^^ii'„(AdSmax)^, m6Si„ax«:l. (36) 

If mbSmax ^ 1, however, we are able to calculate Fq: 



2 — nj \ mb 



2-T1 

■■ + 2 



Fo = 27rA2Ai? ( W„D„ ( ^ ) , mfes^ax » 1. (37) 



where 



1 1 




2' 2 






2 - 


(- 


3(n 





(38) 
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and 



/3„ = 



n + 2 
2n2 



2(n+l) ( ( i 



,'1 1 
vi/( +_ 

. 2 n 



^' ( • ) is the Digamma function. 

If y = exp(A"/(/)") then for Add < 1 so that (n + l)h{Add)/n < 1 where h is given by Eq. (29), we found that: 



dF^{d) 
dA 



n + 1 
nh{Add) 



0(1/ f) 



Since 2{n + l)/{nh) < 1, dF^/dA drops off faster than 1/d for all n. For {n + l)/{nh{Add)) < land mj^ -C d < \ 
the total force between a sphere and a plate is: 



An^A'^R 
n^h{Add)—^d 



1 + 3 +0{l/h^) 



nh{Add) 



(40) 



If ArfSmax «C 1 or mj, / Ad^ 1 then the Fi term is negligible relative to the d-dependant term. If however A^Smax ^ 1 
and nib /Ad 1, f'l = -Pb(smax, i^b) as given by Eqs. (36 and 37). 

For comparison, note that, the total Casimir force between a sphere and a plate is 



TT^R 

360(i3 ■ 



(41) 



V. PREDICTIONS AND CONSTRAINTS 



We now use the results derived above to make specific predictions and derive constraints on theories with either 
V{(j)) = A;^(l + A"/0") or y(0) = A;5cxp( A" /(/)"); Aq = 2.4 x IQ-^eV. The simplest and most natural scenario 
is A « Aq. With these potentials, the energy density of the chameleon field can be identified with dark energy 
cosmologically. If A = 2.4 x 10"^ eV then 

Ad = d/ (82.2 iim) , A^ = 6.92 x 10"^ /xdyne cm'^ = 6.92 x 10"^ mPa. 

For either of the potentials given above one must require n > or n < —4 for a valid chameleon theory to emerge 
[3]. When ^ d <C ti^^, Eq. (22) is actually exact for theories with power-law potentials. In these power-law 

theories, the force per unit area between two parallel plates when m^^ <C c? <C nij^^ is given by Eq. (26). When 
f^bd > 5, dF^/dA is given by Eq. (24). The total chameleon force in the sphere-plate geometry is (when mbd <C 1) 
given by Eq. (34) if n > 2 orn < -4 or by Eq. (35) if < n < 2. 

We begin by considering how Casimir force measurements presently constrain chameleon theories, and then discuss 
the prospects for the detection of chameleon fields by the next generation of such tests in the next section. 

The first attempt to measure the Casimir force between two parallel plates was made in 1958 by Sparnaay [15]. 
The data he found contained large systematic errors, due mostly to the determination of d, and so was only said to 
"not contradict Casimir's theoretical prediction". At the largest separations probed (^ 2/im) this experiment was 
sensitive to pressures between the two plates of 0.1 mPa (=1 mdynecm"^), however the inaccuracy in determining d 
was generally ±0.12 /xm. This measurement was conducted in a vacuum with pressure 10"^ torr. 

The Casimir force between two parallel plates was successfully measured by Bressi et. al. [16]. They measured the 
Casimir force between 0.5 — 3.0 y^m to an average precision of 15%. This corresponds to a sensitivity of approximately 
1 mPa. A vacuum pressure with 10"'^ torr was used in this experiment. 

The most accurate measurements of the Casimir force between two parallel plates have, however, been made by 
measuring the gradient of the force between a sphere and a plate. Dynamical measurements of the force between a 
sphere and a plate would detect not F^°*'{d) but dF^°^/dd and hence, by Eq. (33), dF^/dA. 

We define = dF^/dA to be the chameleonic pressure between two parallel plates. The Casimir pressure between 
two such plates is similarly defined to be Pc = tt"^ /240d'^. Thermal corrections to the Casimir force [21] are sub-leading 
order at the separations that have been probed thus far, and so we do not consider them at this point. To date, the 
most accurate measurements of Pc over separations d ~ 0.16 /um — 1.2 /zm have been made by Decca et al. in a series 
of three experiments taking place between 2003 and 2007 [17-19]. We define P to be the total measured pressure 
between two parallel plates. Using their most recent experiment, described in Ref. [19], Decca et al. found the 
following 95% confidence intervals on ZAP = P - P^: aX d = 162nm, |AP| < 21.2mPa, at d = 400nm, |AP| < 0.69mPa 
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FIG. 2: The solid lines in Figure (a) show the predicted chameleonic pressure between two parallel plates for V = 
Aq (1 + A"/<^"), n > and A = Ao = 2.4 x 10~''eV. The dotted lines show the current experimental constraints on any 
such pressure. Sparnaay, Padova, IndianaOS and IndianaO? refer to Refs. [15], [16], [17] and [19] respectively. The predictions 
shown in Figure (a) only apply when the test masses have thin-shells and for ^ d <C "^j^^; mc is the chameleon mass 

inside the test masses and rrib is the chameleon mass in the background. The white region in Figure (b) shows the values of the 
chameleon to matter coupling, M, for which the predictions shown in Figure (a) are applicable to the most recent experiment 
conducted by Decca et al, labeled Indiana07 in Figure (a). 



and at d = 746nm, |AP| < 0.35mPa. In the first experiment [17], measurements were also made for larger separations. 
For 450nm < z < 1200nm, they found | AP| < 0.54mPa. At d = 162nm, 400nm and 746nm the results of Decca et al. 
represent a detection of the Casimir force to an accuracy of 0.19%, 0.9% and 9.0% respectively. A vacuum pressure 
of 10~^torr was used in making all of these measurements. 

In FIGs 2a and 3a we plot vs. d for <^ d m^^ as solid lines for representative values of n: n — 1/2, 1, 4, 10 
in the former plot, and n = —4, —6, —8, —10 in the latter for V{(f>) — Aq(1 + A"/(/)"). In all these plots we have taken 
A = Aq = 2.4 X 10~^ eV. The dotted lines show the experimental limits on |AP(2:)| and the labels Sparnaay, Padova, 
IndianaOS and Indiana07 refer to Refs. [15], [16], [17] and [19] respectively. For Padova we have taken the upper 
bound on \AP\ to be ImPa. 

It is very clear from FIG. 2a that the magnitude of the chameleonic pressure, P^, predicted by theories with n > 
and A = Ao = 2.4 x 10^"^ eV currently lies well below the experimental limits. The predicted Pc{d) is everywhere 
at least 2 orders of magnitude smaller than the current experiment bounds. The story is very different for n < —4 
theories with the same value of A. FIG. 3a clearly shows that the n = —4 and —6 theories are strongly ruled out 
by the latest 95% confidence limits found by Decca et al. [19] (labeled Indiana07 on the plot). Indeed the n — —4 
theories is even ruled out by the 1958 measurements made by Sparnaay [15]. The value of predicted by the n = —8 
theory is close to edge of what is currently allowed. For n > 0, the larger n is, the steeper the drop-off in P,f,{d) with 
d and, as a result, the larger P^[d) is at separations < A^^ w 82 ^m. Very shallow potentials (0 < n < 2) predict 
the smallest P,p{d). This is disappointing, as the shallower the potential is, the larger the difference (A(/)) between 
the value of (p here on Earth and in the cosmological background. Both variations in the traditional 'constants' of 
Nature and the magnitude of the chameleon force between distant thin-shelled bodies grow with A(/) [3] . The larger 
A(f> is, then, the more scope there is for the presence of a chameleon field to produce non-negligible and potentially 
detectable alterations to the standard cosmological model. 

The values of P^ plotted in FIGs. 2a and 3a are accurate provided d <^ and the test masses have thin- 

shells. It is clear that the strongest constraint comes from the 2007 experiment of Decca et al.. For these constraints 
to actually be comparable with the plotted values of Pc, it must therefore be the case that <C d ^ m^^ for 
d ~ 0.2 — 0.8 fim. We must also require that the test bodies used in this experiment have thin-shells. If Al is too 
small (/3 too large) then rribd ^ 1 for the above range of separations. The chameleonic force between the plates would 
then be exponentially suppressed (by a factor « exp(— TOf,d)) and as such would be negligible. If M is too large (/3 
too small) then the test bodies will either lose their thin-shells or nicd > 0(1). In the absence of a thin-shell, the 
chameleon force has a Yukawa form with mass mf, and we would certainly have nibd 1. Casimir force experiments 
such as those conducted by Decca et al. are, however, only sensitive to Yukawa forces for which rTifjd ^ 0{1). If the 
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Parallel Plate Constraints on n < -4 theories Applicability of Predictions shown in (a) 




Separation (nm) Slope of Potential (n) 

FIG. 3: The solid lines in Figure (a) show the predicted chameleonic pressure between two parallel plates for V — A^ + A"*"*""/;/!)", 
n < —4 and A = Ao = 2.4 x 10""^ eV. The dotted lines show the current experimental constraints on any such pressure. Sparnaay, 
Padova, IndianaOS and Indiana07 refer to Refs. [15], [16], [17] and [19] respectively. The predictions shown in Figure (a) only 
apply when the test masses have thin-shells and for <C d <C m'^ ; rric is the chameleon mass inside the test masses and 
nib is the chameleon mass in the background. The white region in Figure (b) shows the values of the chameleon to matter 
coupling, M, for which the predictions shown in Figure (a) are applicable to the most recent experiment conducted by Decca 
et al, labeled IndianaO? in Figure (a). 



test-masses do not have thin-shells then, Casimir force experiments cannot be used to constrain chameleon theories. 
If rricd > C(l), then « const and so, once again, Casimir force experiments would be unable to constrain it. FIGs. 
2b and 3b shows the values of M = Mp\/ (3 for which the values of shown in FIGs. 2a and 3a can be compared 
with the IndianaO? constraints. 

Parallel Plate Constraints on V = a" exp(A"/(^") theories; A = 2.4 x 1 o^^eV 
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FIG. 4: The solid lines show the predicted chameleonic pressure between two parallel plates for V = AoexpA"/(^" and 
A = Ao = 2.4 X 10~^ eV. The dotted lines show the current experimental constraints on any such pressure. Sparnaay, Padova, 
IndianaOS and Indiana07 refer to Refs. [15], [16], [17] and [19] respectively. The predictions shown above only apply when the 



test masses have thin-shells, 
in the background. 



is the chameleon mass inside the test masses and mj, is the chameleon mass 



The predictions and constraints found above apply to theories with V{(j)) w Ag(l + A"/0"). li V = Ag(7(A"/0"), 
where g{y) has a Taylor expansion about y = and g'{0) = 1, then the above predictions would apply for (A/0)" ^ 1 
i.e. <C A. We expect A = Aq « 2.4 x 10"'^ eV so that the chameleon field is responsible for the late time 
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acceleration of the universe. The above predictions would then only apply if d ^ « 82 ^m, which is not the 
case. Generally speaking, if the potential is very steep, = i^^/yl cx 1/d^ which is a stronger d dependence than 
that exhibited by theories with V = Ao(l + A"/0") and n > 0, but a weaker d-dependence than that predicted by 
power-law theories with n < —4. 

For concreteness, we consider the parallel plate predictions and constraints for a theory with V ~ AgexpA"/^" 
and again A = Ao = 2.4 x 10^"^ eV. For d ^ 82 /im we find if ti > the chameleonic pressure, Pct,{d), predicted by 
such theories is larger than that predicted by theories with V = Aq(1 + A"/</)"); if n < —4, the opposite is true. FIG. 
4 shows the predictions for made by a theory with = Ag exp A"/(^". It is clear from these plots that currently no 
such models are ruled out, although the predicted chameleonic pressure for n = 1/2 is an order of magnitude larger 
than it would if V{(t)) were exactly A;^(l + A"/0"). 

The experiments performed by Decca et al. employed the sphere-plate geometry but measured the Casimir force 
dynamically. Dynamical experiments such as these directly measure not forces but the rate of change of forces with 
separation i.e. dF/dd. Since dF/dd (x dF/dA, what is actually measured is equivalent to the force per unit area 
between two parallel plates. Other experiments that use the sphere-plate geometry have made static, rather than 
dynamical, measurements of the force between the two bodies. Static measurements allow one to place limits on the 
force itself rather than its gradient. It should be noted that when such measurements are performed, it is necessary to 
calibrate the experiment so as to eliminate any electro-static forces. This has the effect that at some large separation, 
rfcai say, the force between the sphere and the plate is defined to be zero. If one has an expression for the force F{d) 
then, what would actually be measured by these experiments is AF{d) — F{d) — -F(c?cai)- These experiments are 
therefore insensitive to forces that are virtually constant for d < deal- 

In 1997, Lamoreaux measured the force between a spherical lens with radius (12.5 ±0.3) cm and a 2.54 cm diameter, 
0.5 cm thick optical flat [10]. Measurements of the Casimir force where made for separations in the 0.6 to 6 /xm range 
and the Casimir force was measured to an overall accuracy of 15%. The calibration of the system was performed at 
a separation of about 10 /im. At the largest separations (d > l/xm), the results of this experiment place an upper 
bound on the magnitude of any residual force, which includes any thermal corrections to the Casimir force, of about 
30 pN. The experiment was conducted in a vacuum with pressure 10-''torr. 




FIG. 5: The solid lines in Figure (a) shows (_F^°*(d) — F^°^ (deal)) / R where F^°^ is the predicted chameleonic force between a 
sphere and a plate, for V = Aq{1 + A" /(j)"), n > and A = Ao = 2.4 x 10"'^ eV. R is the radius of the sphere and we have taken 
deal = 10 /im. The dotted line shows the current best experimental constraint on any such force pressure, which comes from 
Ref. [10]. The predictions shown in Figure (a) only apply when the test masses have thin-shells, rricd 2> 1 and mtdcai ^ 1; rric 
is the chameleon mass inside the test masses and nib is the chameleon mass in the background. The white region in Figure (b) 
shows the values of the chameleon to matter coupling, M, for which the predictions shown in Figure (a) are applicable to the 
1997 Casimir force measurement performed by Lamoreaux [10]. 



A similar measurement was made in 1998 by Mohideen & Roy [13]. In this experiment a relatively small polystyrene 
sphere was used with diameter 196 /xm, and the system was calibrated at a separation of about 900 nm. Mohideen 
& Roy were able to measure the Casimir force to precision of 1% at the smallest separation of about 100 nm. They 
found an RMS derivation between experiment and theory of 1.4 pN [13]. The error bars on the measurements at 
individual separations were however larger, being about ±7 pN. This experiment was performed in a vacuum with 
pressure SOmTorr. 
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In additional to dynamical force measurements, in their 2003 experiment Decca et al. also made a static measure- 
ment of the Casimir force [17]. The experiment was calibrated at a separation of 3 /xm. The direct force measurements 
limit |Ai^| < 0.5 pN for 400 nm < d < 1200 nm. In this experiment the sphere had radius (296±2) fim and the vacuum 
pressure was 10~^ torr. 

We found in Section IV B that when the sphere and the plate have thin-shells and m^^d <C m^^ then if deal <C tIj^^: 



where — Aq/A. If, however, deal ^ ttIj^ ^ then 



2'nRKlKKn 



(Add)- — 



Dr. 



rrib 



It should be noted the terms that for 



c^cai 3> d, the terms that depend on mi, and deal a-re only important for 



< 71 < 2 theories. For all n, the larger R is (for fixed d), the larger AF^°*(c?). R is largest (by several orders of 
magnitude) in the 1997 Lamoreaux experiment [10]. The relatively large dimensions 0{cm)) of the test masses 
used in this experiment (see Ref. [10]) ensure that they have thin-shells for a larger range of M than do the test masses 
used in Refs. [13] and [17]. Currently then, the best constraints on chameleon theories from static measurements of 
the Casimir force using the sphere-plate geometry are provided by the Lamoreaux experiment. 



Sphere-Plate Constraints on n >0 theories 
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FIG. 6: The solid lines in Figure (a) shows (F^°'(d) — _F*°*(dcai))/i? where F^°^ is the predicted chameleonic force between 
a sphere and a plate, for V = Ao(l -|- A"/(^"), n < —4 and A = Ao = 2.4 x 10~^ eV. R is the radius of the sphere and we 
have talcen deal = 10 fim. The dotted line shows the current best experimental constraint on any such force pressure, which 
comes from Ref. [10]. The predictions shown in Figure (a) only apply when the test masses have thin-shells, rricd ^ 1 and 
rribdcai ^ 1; fnc is the chameleon mass inside the test masses and mb is the chameleon mass in the background. The white 
region in Figure (b) shows the values of the chameleon to matter couphng, M, for which the predictions shown in Figure (a) 
are applicable to the 1997 Casimir force measurement performed by Lamoreaux [10]. 



In FIGs. 5a and 6a we plot, as solid lines, the predicted values of AF^°'(d)/_R with deal = 10 /im, such as it is 
in Lamoreaux experiment [10]. In making these predictions, we have taken V{4)) — Aq(1 + A"/(/)") and A = Aq = 
2.4 X 10~^eV. These predictions are accurate provided < d < deal ^ ^n^^ and the test masses have thin-shells. 
If either of these conditions did not hold, then AF^°* would be much smaller. The dotted line in each plot is the upper 
bound placed on AF^"*/R by the 1997 experiment of Lamoreaux [10]. FIG. 5a shows the predictions for n > 0. We 
see that if n > 0, the predicted values of AF^°'^/R are several orders of magnitude smaller than the current experiment 
upper bound. Predictions for n < —4 are shown in FIG. 6a. It is clear to see that theories with this potential and 
A = 2.4 X lO^'^eV and n = —4 and —6 are strongly ruled out by the Casimir force measurements made by Lamoreaux 
[10]. This picture is remarkably similar to that found by comparing measurements of the force between two parallel 
plates with the predictions of chameleon theories; there too theories with n = —4 and —6 were ruled out. 



19 




FIG. 7: The solid lines show {F^°^{d) — F^°^ (dcsd)) / R where -F^°* is the predicted chameleonic force between a sphere and a 
plate, for V — AQexpA"/(^" and A = Ao = 2.4 x 10~^ eV. R is the radius of the sphere and we have taken deal ~ 10 ^m. 
The dotted line shows the current best experimental constraint on any such force pressure, which comes from Ref. [10]. The 
predictions are valid for ^ d < deal ^ "^(T^i ^ ^^e chameleon mass inside the test masses and nit is the chameleon 
mass in the background. 



For the predictions shown in FIGs. 5a and 6a to apply to the 1997 experiment performed by Lamoreaux, we must 
require that the test masses have thin-shells, rricd 3> 1 and that TOfcdcai <C 1. If either of these conditions fail then 
the magnitude of the force due to the presence of the chameleon field would be much smaller than the predictions 
shown. The white region in FIGs. 5b and 6b indicates the values of the chameleon to matter coupling, M, where 
these conditions are predicted to hold for an experiment such as Lamoreaux's [10]. FIG. 7 shows how the predictions 
for AF^°* = (d) - (deal)) with V ^ exp A'' and A = Aq = 2.4 x lO'^ eV compared to the experimental 
data. It is clear that currently no such theories are ruled out. 

We conclude this section by presenting the current combined constraints from Casimir force experiments on 
Chameleon theories with V{cf)) = Ag(l + A^/c^") and matter coupling M; Aq = 2.4 x 10"'^ eV. The current combined 
constraints on M and A for theories with n = —8, n = —4, n = 1 and n — 4 are show in FIG. (8). As we noted 
above chameleon theories with n = —4, A w Aq and M < 5 x lO^'^ GeV are currently ruled out. The n = —4 
corresponds to a theory where V{(t)) = Aq + Xcj)"^ where A — (Aq/A)'*. This potential first considered in the context of 
chameleon theories in Ref. [20]. The strongest constraint on A from Casimir force measurements is A < 2 x 10~'*eV 
i.e. A > lO"' for 10** GeV < M < 10*^ GeV. In the n = -8, n = 1 and n = 4 plots we can see that A « Aq is not 
currently excluded (indeed this the case for all theories with n < —8 or n > 0). In general, the steeper the drop-off 
of F^/A with d is, the stronger the constraints on A are. For n > 0, the steepness of F^/A increases as n — > oo; for 
n < —4 the steepness decreases as n ^ — oo. Currently the strongest constraints on A are: for n = 1, A < 40 eV for 
2 X 10^ GeV < M < 3 X 10*^ GeV; for n = 4, A < 0.1 eV for 2 x 10^ GeV < M < 2 x lO" GeV and for n = -8, 
A < 3 X 10-3 eV for 2 x 10^ GeV < M < 7 x 10^^ Q^y poj. potentials like V = A"^ exp(A/0)", which are steeper than 
any power-law, we generally have A < 10 — IOOAq for 10^ S ^ ^ 10^^ GeV. Although Casimir force measurements 
are generally unable to see chameleon theories with A « Aq at present, they do generally require A be no more than 
a few orders of magnitude larger than Aq. Unless we are prepared to allow for two small but unrelated energy scales 
in the theory, it therefore, if chameleon fields do exist, seems natural that we should expect A « Aq. 



VI. PROSPECTS FOR FUTURE EXPERIMENTS 



Measurements of the Casimir force have so far been unable to rule out or detect chameleon theories where the energy 
scale of the potential is approximately the dark energy scale of 2.4 x 10"^ eV. The exception are the n = —4 and 
n = —6 chameleon theories with Ratra-Peebles potentials. In this section we consider the prospects for future Casimir 
measurements detecting chameleon field theories, in particular those with A — 2.4 x 10^^ eV. We shall identify two 
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Current Constraints from Casimir Measurements on Chameleon Theories Current Constraints from Casimir Measurements on Chameleon Theories 
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FIG. 8: Current constraints from Casimir force measurements on chameleon theories with V{(fi) ~ Ao(l + A"/(/i"); Aq = 
2.4 X 10"^ eV. The cases n = —8, n — —4, n — 1 and n = 4 are shown above. Presently A ~ Ao is only ruled out for theories 
with n — —4 and n = —6 (for which a plot is not shown). See text for further discussion. 



proposed experiments that could potentially make such a detection. 

We define e{d) to be the ratio of the chameleonic, F^, and Casimir, Fcas forces between two parallel plates: 

eid) = 

If the chameleonic force is to be detected or ruled out, one would have to measure the Casimir force at a separation 
d to an accuracy of 100e(d)%. If ¥{(/)) = A^G{{ A /(/))"), for some G such that G(l) = 1 and G"(l), G"(l)/2 ~ 0(1), 
then we would find Ftp/ A « A(f,AQH{Kdd) where h.d = ^0/^ ^'^"^ some function related to G; H{1) = 1 and 
~ 0(1). It follows that 

,,, 24OA0A4 n4„.A n 

For ^ C'(l), 2AQA^/tt'^ ^ 0(20). Thus in the simplest and most natural situation where there is only one energy 
scale associated with the potential i.e. A « Aq, we expect £(Ag ^) ~ 0(20). If G{Add) oc d~^ and A « Aq, then 
£{d) ~ 0(1) when d - 0{Aq'^/3). If we take A « Aq = 2.4 x IQ-^eV then A^ ^ 82 ^m. We would then expect 
e{d) ^ 0(1) when d « 30 /xm. At d « 10 /im, the ratio of the chameleon to the Casimir force would then be a 
few percent. This provides us with a rough estimate for the sensitivity required to detect chameleon theories with 
A « Aq « 2.4 X 10^'^ eV. A more precise requirement can, of course, be given when V{(l)) is specified. 
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Relative Strengths of Forces between two Parallel Plates 
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FIG. 9: Relative strengths of the chameleonic and Casimir forces between two parallel plates at zero temperature. We have 
taken A = 2.4 x 10"^ eV. Figure (a) is for V ^ A^^ + A'^+"/(j>" and figure (b) for V = A"^ exp(A"/<^") with A = 2.4 x 10"^ eV. 
For most values of n we see that e{d) = 1 for d ~ 30 — 40 ^m. At d ~ 10 /im, e ~ 0.01 — 0.1 in most cases. 



If V" = (1 + A"/(/)") then for rribd < 1: 



'2d()K A** 2{n+4.) 
e{d) = epow(d; n) = (Add)^^ . (42) 



If y = A(^exp A'7(A" then for nibd < 1 and Aod < 1 so that (n + l)/nh{Ad) < 1: 

/ , N 480A4 n + 1 

e{d) « ecxp(d; n) = (Arfd) 
n2A5/i(Arfcf)T~ 



nh{Add) 



(43) 



If F = AqCxp A"/0" and d A^ ^ we have e — Epow For both of these examples F^/A grows more slowly than the 
Casimir force as d — > 0. We plot e{d) against d for both of these potentials in FIG. 9. In both figures we have taken 
A = Ao = 2.4 X 10"3 eV. For most values of n we see that e{d) = 100% for d « 30 - 40 ^m, and e{d) « 1% - 10% when 
d w 10 /X. In principle at least, all of these chameleon theories could be detected if one were able to unambiguously 
measures forces that were 1% the size of zero-point Casimir force at cZ « 10 fim. This corresponds to a sensitivity to 
pressure, P = F/A, between the two plates of 0.13pNcm~^ = l.SnPa. As we discuss further below, such a precision 
is well within the reach of the next generation of experiments. However, in order to actually detect the chameleonic 
force one must not only be able to reach this sensitivity, but also be able to calculate and control all background 
non-chameleonic forces to the same precision. Whilst it is possible to do this for the zero-point (i.e. zero temperature) 
Casimir force and any electrostatic forces, it becomes a problem at room temperature and at separations of about 
10 fim. Under these conditions the thermal contribution to the Casimir force is expected to dominate over the zero- 
point force. If one were able to calculate the thermal contribution to the Casimir force accurately then this would not 
be a major problem. However, there are currently two main approaches to calculating the thermal Casimir force and 
there is a great amount of debate and disagreement as to which is correct [21]. The two models are often refereed to 
as the Drude model [22] and the plasma model [23] . For an excellent review of the current status of this controversy 
see Ref. [21]. 

For separations larger than S/im, the thermal Casimir force between two perfectly reflecting mirrors is: 

-^thermal _ C{3)kBT 

A " And^ ' 

where is Boltzmann's constant. The competing more realistic models predict FthGrmai/2 < ^thermal ^ ^thermal- At 
T = 300K and at d = 10 ^m we then have 1.5 < -Fthormai/^cas ^ 3.0. At 10 /^m , the magnitude of the chameleonic 
force for A = Aq ~ 2.4 x 10~^eV is generally about 1% of the total Casimir force (including thermal correction). 

In FIGs. 10 and 11 we indicate how F^/{Fcas + -^thermal) and f^/f thermal depend on T and d. 

To detect, or rule out, the presence of chameleon fields with A w Ag « 2.4 x 10""^ eV at a separation of 0(10) /im, 
one first has to calculate both the zero-point Casimir force to an accuracy of better than 1% and then either to be 
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Dependence of Force Ratio for two Parallel Plates on Temperature °' ^"'^ '^thermal d - 10n 
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FIG. 10: Relative strengths of the chameleonic and the Casimir forces between two parallel plates at different temperatures. 
We have taken V = Ao(l + A")/^" and A = Ao = 2.4 x 10"^ eV. Figure (a) shows how F^/{F^^s + -Flhermai) at d = 10 /im 
depends on T, and figure (b) shows the temperature dependence of -F^/^thermai. If n = 1/2 then we see that the chameleonic 
force only dominates over the thermal contribution to the Casimir force for T < 2 K. Figure (c) shows how / {Fcas + -^thermal) 
at 300 K depends on separation, d. We see that at 300 K and with d « 10 /im, e ~ 0.005 — 0.1. 



able to do the same for the thermal Casimir force. Otherwise perform the experiment at sufficiently low temperature 
so that Fthermai ^ F^. From FIGs. 10 and 11 we see that to realize this latter option for most theories considered 
here, the experiment would have to be run at T ~ 0{1)K or cooler. 

The Gravitation group at the University of Birmingham [24] have constructed a super-conducting torsion balance 
that they intend to use to measure, amongst other things, the Casimir force at 4.2 K [24]. Precise details of the 
separations they will probe and the precision which they expect to achieve have not yet been announced. However, 
even at ci « 10 /Ltm, it is clear from Figures 10b and lib that even with T = 4.2 K the chameleonic force between two 
parallel plates is still generally only 0.3 — 3 times the size of the thermal contribution to the Casimir force. Even at 
this low temperature, one would still have to be able to calculate the thermal Casimir force to greater accuracy than 
is currently possible due to the controversy over the different models. 

The thermal Casimir force drops off as l/cP as d ^ oo. Of the chameleon theories considered here only those 
with n = —4 and ri = —6 exhibit a faster drop off for d ^n^^ and with a Ratra-Peebles power-law potential, 
these theories are already ruled out for A w Aq = 2.4 x lO^^eV. If one wished to avoid having accurately model the 
thermal Casimir force, one could take advantage of the slow drop-off of the chameleonic force, and run a Casimir force 
experiment at separations where the chameleonic force is predicted to dominate over the total Casimir force. For 
A « Aq « 2.4 X 10"'^ eV, this would generally involve running the experiment at separations larger than 30 — 100 ^m. 
Although this range of separations has already been probed by tests of gravity such as the Eot-Wash experiment 
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FIG. 11: Relative strengths of the chameleonic and the Casimir forces between two parallel plates at different temperatures. 
We have taken V = Aq exp(A"/(/)' and A = Ao = 2.4 x 10~^eV. Figure (a) shows how f^/(Fcas + -Fthermai) at d = 10 ^m 
depends on T, and figure (b) shows the temperature dependence of -F^/-Fthcrmai • If n = 1/2, the chameleonic force dominates 
over the thermal contribution to the Casimir force for T < 5K. Figure (c) shows how F^/{Fcas + -Fthermai) at 300 K depends 
on separation, d. At 300 K and with d 10 //m, e ^ 0.002 - 0.02. 



[14], the physical electrostatic shield that is employed in such tests would block any signal due to strongly coupled 
chameleon fields [3] . 

Importantly, the way in which electrostatic forces are controlled in Casimir force measurements, i.e. without a 
physical shield, does not shield the chameleonic force [3]. Hence any such experiment could probe relatively large 
separations and could use large test masses (i.e. with lengths scales of a few centimeters or more) so as to magnify 
any new forces. Given the size of the test masses, and the very high accuracy force measurements that would be 
required, it would seem sensible for any such experiment to make use of a torsion balance as was used by Lamoreaux 
in Ref. [10]. Lamoreaux has discussed a number of improvements to his 1997 experiment [25] which would allow it to 
detect chameleon fields with A « Aq w 2.4 x 10^'^ eV. We discuss this further below in section VI B. 

In the longer term, there are a number of planned experiments (not least the new one proposed by Lamoreaux) that 
aim to detect the thermal Casimir force to high accuracy. If experiments can decide which is the best model, then 
one of these new experiments, proposed in Ref. [26] and already under construction at the Institute Laue-Langevin 
(ILL) in Grenoble, would ultimately be able to detect or rule out the presence of chameleonic forces at d w 10 ^vci. 
We discuss this experiment further below. 
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A. Proposed Grenoble Experiment 



In Ref. [26], Lambrecht et al. proposed a new experiment to measure the Casimir force between two parallel 
plates and search for thermal corrections in the 1—10 fim separation range using a high sensitivity torsion balance. 
Previous attempts to make such measurements [15, 16] using this geometry have been limited by the extent to which 
the plates can be kept parallel. The solution that Lambrecht et al. proposed in Ref. [26] was to take advantage 
of the inclinometer developed for a neutron experiment performed at the Institute Laue-Langevin (ILL) in Grenoble 
[28]. This would allow them to limit deviations from parallelism at the 10~^radians level. For comparison the Padova 
experiment of Bressi et al. [16] had estimated deviations from parallelism of 3 x lO^^radians. The surface area of 
the plates that will be used in the new experiment would be of the order of 120 cm^. The force measurements made 
using the torsion balance would have a resolution of about IpN for 5 /im ^ d < 10 /im [26]. The aim of this proposed 
experiment is two-fold. Firstly, it should contribute to the settling of the controversy that surrounds the correct 
approach to calculating the thermal Casimir force. Secondly, with a theoretical model for the thermal Casimir force 
established, it would be able to probe for the presence of new forces. 
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FIG. 12: Shaded area is the region of A-n parameter space that could potentially be detected or ruled out with 95% confidence 
by the new Grenoble experiment proposed by Lambrecht et al. in Ref. [26] . Figure (a) is for V = Ao(l + A")/(/>"- and figure (b) 
for V = Aq exp(A"/(/!>") with Ao = (2.4±0.1) x 10"^ eV. If A = Ao, which would be natural is the chameleon field is responsible 
for the late time acceleration of the universe, then A must lie between the two dotted lines i.e. A = (2.4 ±0.1) x 10"'^ eV. It 
should be noted that for these constraints to apply the chameleon to matter coupling /3 = Mpi/M must be large enough for 
test masses to have thin-shells and small enough for the inverse chameleon mass in the background to be large compared with 
the separations probed i.e. m^"^ ^ lOfim. Both of these requirements will depend to varying degrees on the quality of the 
laboratory vacuum in which the experiment is performed. 



In order for there to be a detectable chameleonic force one must ensure that the laboratory vacuum is good 
enough so that the chameleon mass in the background, mb, is <C hc/lOfj.m. Provided that this is the case and the 
IpN sensitivity can be reached with background and systematic effects controlled at the same level, the Grenoble 
experiment proposed in Ref. [26] would, by making measurements for distances 5^m < d < 10 fim be able to detect, 
or rule out, a large number of chameleon theories. Specifically ii V — AQexp(A"/0") and A = Aq = 2.4 x 10^"^ eV, 
then this new experiment could detect or rule out all theories with 1/8 < |n| < 80 with at least 99.5% confidence. 
Additionally, it would be able to rule or detect all theories with V{(j)) = Ag(l -I- A")/</)", A = Aq = 2.4 x 10"^ eV and 
n < —4 or n > 1/3 with the same confidence. Figure 12 shows 95% confidence limits on the values of n and A that 
this new test should ultimately be able to place on chameleon theories. 

The region of M space that will be accessible to this new experiment will depend on, amongst other things. A, n 
and the quality of the laboratory vacuum and the thickness and composition of the test masses. Firstly, it is required 
that the test-masses have thin-shells. In this experiment the test masses are 15 mm thick glass plates. The thin-shell 
conditions generally hold for all /3 = Mp\/M > 0(1). 

Constant forces are not detected by this experiment. We define nic to be the chameleon mass inside the test 
masses. If nicd <C 1 for 5/j,m < d < lO^m then the force between the plates will be virtually constant, and therefore 
undetectable. Additionally if mbd ^ 1, the chameleonic force will be undetectably small. 
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Potential Constraints from new Grenoble Experiment with P =10 torr 
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Potential Constraints from new Grenoble Experiment with P =10"^ torr 
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FIG. 13: Shaded area is the region of A-M parameter space that could potentiaUy be detected or ruled out with 95% confidence 
by the new Grenoble experiment proposed by Lambrecht et al. in Ref. [26] for two choice of the vacuum pressure. The plots 
shown are forn = 1 with either V = A^(l + A")/<?!)" or F = exp(A"/,^"). The situation for other 0(1) values of n is similar. 
Ao — (2.4 ± 0.1) X 10"'' eV. If Pvac ~ 10~^torr then for both choices of potential, the Grenoble experiment could detect all 
n = 1 theories with A = 2.4 x 10** and lO^'GeV < M < 10^'^GeV. If Pvac is lowered to lO^'torr, then values of M smaller than 
lO'^GeV could be detected. 



Figure 13 shows the region of A/-A parameter space that the new Grenoble experiment could potentially detect for 
n — 1 and two choices for the vacuum pressure, Pvac- We have taken V = Ao(l + A^/cj)") and n ~ 1. The picture is 
very similar for other 0(1) values of n and for V = Aq exp(A"/^"). It is clear that the better the vacuum pressure, 
the smaller the values of M that can be detected. If Pvac = 10~^torr then for both choices of potential, the Grenoble 
experiment could detect all n = 1 theories with A = 2.4 x lO'* and lO'^GeV < ^/ < lO^^GeV. If Pyac is lowered to 
lO^torr, then values of M smaller than lO'^GeV could be detected. 

Since values of M smaller than lO^'GeV are not generally consistent with particle physics [6], a lO^torr vacuum 
would be sufficient to detect or rule out all chameleon theories with V — Aq{1 + A" / cf)) or V — AQexp(A"/0") and 
A = Ao = (2.4 ± 0.1) X lO-^eV and M < 10^^ GeV. 

Provided the laboratory vacuum is of sufficient quality, the Grenoble experiment proposed by Lambrecht et al. 
[26] has the potential to detect or rule out almost all of the strongly coupled chameleon theories with A = Ao = 
(2.4 ± 0.1) X 10~^eV consider here. Before this could be done however, the controversy that surrounds the correct 
method for calculating the thermal Casimir force would have to be settled, and the total Casimir force modeled 
theoretically to an accuracy of better than 1% at d ~ 5 — 10 /im. 
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B. New Lamoreaux Experiment 

From Figures 10c and 11c it is clear that the chameleonic force predicted by n > theories with either V = 
Aq(1 + A"/0") or V = Agexp(A"/0") and A w Aq = 2.4 x 10~^eV dominates over the expected thermal correction 
when d « 70 iim — 300 iim. At such large separations the chameleonic force between two parallel plates is more than 
100 times the size of the zero-point Casimir force. If one wished to detect or rule out chameleonic forces but to avoid 
the controversy that surrounds the thermal Casimir force, one could therefore search for new forces at separations 
where the chameleonic force dominates over the total Casimir force. 



Potential Constraints from New Lamoreaux Experiment Predictions for New Experiment Proposed by Lamoreaux 
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FIG. 14: Shaded area is the region of M-A parameter space that could potentially be detected or ruled out the new experiment 
proposed by Lamoreaux [25, 27] i.e. AF^°'(33Atm) = _F^°'(33Atm) - -F^°'(dcai = 1cm) > 0.1 pN. The plot on the left is for 
V = Ao(l + A")/(^", Ao = 2.4 X 10"^ eV and n = 1, and the one on the right is for V ^ At exp(A"/0"), Ao = 2.4 x 10"^ eV 
and n — 1. The graphs for other values of n close to 1 are similar. 



Lamoreaux has proposed a number of improvements to his 1997 torsion balance experiment [25, 27]. The new 
experiment would be stable for separations as large as 1 cm [25] . The improved experiment would make use of one 
flat plate and a curved one with radius of curvature 3 — 4 cm[25]. The pressure of the vacuum is intended to be 
5 X 10^^ torr [25] . Ultimately, this new experiment should be able to detect changes in forces that are as small as about 
O.lpN [25]. We assume that the experiment is electrostatically calibrated at deal ~ 1cm, and that the chameleonic 
force can be distinguished from any residual electrostatic forces (which would behave as 1/d and 

We define AF^°\d) = F*°*(d) - F*°*(l cm). The total Casimir force is less than 0.05 pN for d > 33 fim. Therefore, 
even without a precise knowledge of the thermal corrections to the Casimir force, the new test proposed by Lamoreaux 
should be able to detect the chameleonic force provided AF^°'(d = 33 /xm) > O.lpN. 

Figure 14 shows the region of the M — A parameter space of n = 1 chameleon theories that the new Lamoreaux 
experiment will be able to detect or rule out i.e. AF^ i.e. the values of M and A for which AF^{d = 33 /Ltm) > 0.1 pN. 
We have plotted for both V = Aq{1 + A"/(^") and V = Aq exp(A"/(^"). We note that for n = 1, this new test would 
be able to detect 10^ GeV < M < lO^^ GeV for A = Ao. Furthermore, it would be sensitive to values of A as small as 
7 X 10~^eV. With either potential, for n > 0, the chameleonic force between a sphere and a plate drops more slowly 
than 1/d. 

For n — 1 and A « Aq = 2.4 x lO^^eV, larger values of M (i.e. smaller chameleon to matter couplings) should be 
detectable by the new Lamoreaux experiment but not the new Grenoble test [26]. However the new Grenoble test 
is sensitive to smaller values of A. For the Grenoble test to detect n = 1 chameleon fields, we found that one would 
first have to be able to model the thermal Casimir force to an accuracy of better than 1%. In the new Lamoreaux 
experiment, however, the chameleonic force dominates over the total Casimir force when d = 33 /im and n — I, and 
so useful constraints could be derived without a detailed knowledge of the thermal contributions to the Casimir force. 

In the Grenoble experiment, a relatively large range of n could ultimately be detected. If one wishes to avoid having 
to deal with thermal corrections, however, this is not the case for the Lamoreaux experiment. In Figure 15 we show 
the predicted contours, in n — M parameter space, of AF^^Aq"^ « 82 /xm) with A = Aq. We have only plotted AF^ 
contours for values of n — M for which the test masses have thin-shells. In the absence of a thin-shell, the chameleon 
field behaves simply as a Yukawa field with mass m?,. The 0.1 pN contour represents the limit of detectability for 
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FIG. 15: Shaded areas show the contours of the predicted values A_F|°' (82 /^m) = F^°*{82fim) — _F0°*(dcai ~ 1cm) in the 
new experiment proposed by Lamoreaux [25, 27]. We have taken V{(j>) = Ao(l + A"/(/!>") and A = Aq = 2.4 x 10~^ eV. The 
experiment should be sensitive to AF^°^ > 0.1 pN. At this separation the total Casimir force is < 7.6 x 10^'^ pN, and so, with 
improved sensitivity, it should ultimately be able to distinguish chameleonic forces as 0.01 pN from the Casimir background 
without an incontravertible model for the thermal contribution to the Casimir force. 



the new experiment proposed by Lamoreaux. At a separation d — Aq « 82 /im, the total Casimir force (including 
thermal correction) between the sphere and the plate is less than 7.6 x 10^"^ pN. With the sensitivity of 0.1 pN, the 
Lamoreaux test would be unable, without an undisputed model for the thermal Casimir force, to detect n < 1.6. If 
the sensitivity of the torsion balance could be improved by an order of magnitude, the Lamoreaux would be capable 
of detecting or ruling out chameleon theories with n < 6. Note that, at this separation Ai^^°*, for n < — 8 theories 
would still be too small to detect. The Lamoreaux experiment is most sensitive to chameleon theories with n ~ C(l); 
this is because d-F'^°'/dc? has a slow (but not too slow) drop-off as c? — > oo in these theories. 

VII. CONCLUSIONS 

We have investigated the possibility of using experiments which measure the Casimir force to constrain theories 
where a scalar field is coupled to matter and where there is a chameleon mechanism. Our primary aim was to extract 
the bounds that Casimir force tests currently place on chameleon theories and to make predictions for what near 
future Casimir experiments will be able to detect. 

Chameleon theories are particular interesting because, for certain choices of potential, they can be agent responsible 
for the late-time acceleration. For chameleonic dark energy, one must generally require a potential of the form 
V = A^F{(I)/A) where Aq = (2.4 ± 0.1) x lO'^eV, F = 1 cosmologically today and F'(l),i^"(l) - 0(1) sets the 
scale of A. To the best of our knowledge, chameleon models of dark energy almost always feature at least one scale 
energy scale, Aq, and so do not in themselves alleviate the fine-tuning problems associated with dark energy. If these 
models are to be seen as in some sense 'natural' then one would probably have to postulate the existence of new 
physics at energy scales of the order of Aq. It is clear from this and previous works that chameleon theories with 
A « Aq = (2.4 ±0.1) X 10"'^ eV could have remained undetected thus far are not generally ruled out. This, combined 
with the lack of a generally accepted solution to the naturalness problem of the dark energy scale, means that we 
should not discount the possibility that there really is some new physics associated with the meV scale. Albeit without 
any obvious connection to chameleon fields, this possibility has also be raised in the context of the super-symmetric 
large energy dimensions (SLED) proposal [8]. What makes the SLED proposal and chameleonic dark energy model 
so interesting is that, unlike so many other explanations of dark energy, they make predictions that are eminently 
testable and falsifiable by near future laboratory experiments. If all we ever learn about dark energy comes from 
astronomical observations then cosmic variance alone means it may be difficult to ever fully understand its behaviour 
and its origins. If dark energy can be detected under the controlled conditions of laboratory experiments, however, 
then the prospects for understanding it could potentially be much better. 

In this work, we have aimed to remain as general as possible in our treatment of chameleon theories although when 
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specific predictions have been required we have taken the potential to be of the form of either V{(j)) = Fi = Ag ^1 + ^ j 

or V{(j)) = V2 — Aq exp(A"/(/)"). For different vahies of n and A these potentials were shown to result in a wide range 
of different predictions for Casimir force experiments. Additionally, one can generally think of a more general potential 
as behaving locally as one of those two potentials for some n, A and Aq. These choices for the form of V are therefore 
very handy for understanding the extent to which Casimir force experiments constrain general chameleon theories. 

We found that the magnitude of the chameleonic force that would be detectable by Casimir force measurement, 
generally constrains A < 0(100) Ao = (2.4 ±0.1) x IQ-^eV. Although for very shallow potential larger values of A arc 
still allowed, we still found that generally both A and Aq must be small in experimentally viable chameleonic dark 
energy models. If one small energy scale in the Universe seems undesirable then two unrelated small energy scales is 
at least doubly so. It therefore seems natural that A and Ao be related i.e. A w Aq. We foimd that the chameleonic 
force predicted by theories with A = Aq = 2.4 x 10^'^ eV and V ^ Vi and n > 0, or n < — 10 or ^ = V2 with any n, 
lies well below the present days experimental limits set by measurements of the Casimir force. However, this is not 
the case for theories with V = Vi and n = —4 or n = —6 (and A — Aq) which are strongly ruled out by the latest 
95% confidence limits found by Decca et al. [19]. Additionally, the prediction for V = Vi, n = —8 theory lies close to 
the edge of what is currently allowed. In all allowed theories, the chameleonic force between two parallel plates with 
separation d increases more slowly than as d ^ 0. 

Casimir force experiments are generally more sensitive to largo gradients in forces than they are to large but 
slowly varying forces. As a result, there are only able to place constraints on Yukawa fields with mass, m^, if they 
probe separations, d, such that m^d ~ 0{1). If the test masses used in these tests do not have thin-shells, then the 
chameleon field behaves as a Yukawa field with mass nib, where nib is the chameleon mass in the background. If we 
take TOc 3> rub to be the chameleon mass deep inside the thin-shelled tost masses, then for a largo range of separations 
^ d < m^^ we found that the chameleon force increased as d ^ like some inverse power of d; crucially 
this is a much steeper variation of the force with d than is predicted for a Yukawa field. For nicd < 1, however, 
the chameleonic force between thin-shelled bodies depends only very weakly on d which makes it difhcult to detect. 
If rubd 2> 1 then the chameleonic force is exponentially attenuated. Casimir force measurements therefore provide 
the strongest constraints on chameleon theories for which m^^ ^ d ^b^ which limits the range of chameleon to 
matter couplings that they can detect. Specifically, if A = Aq and the matter coupling has gravitational strength 
i.e. M ~ 0{Mpi) then it is generally the case that nicd < 1 and / or the test masses do not have thin-shells. As a 
result, wo saw that Casimir force tests are best siiitod to soarciiing for and limiting the properties of strongly coupled 
chameleon fields i.e. Mp\/M ^ 1. This is useful, as it is precisely these strongly coupled theories that are the most 
difficult to constrain [3]. If M > O(Mpi) then excellent constrains on chameleon fields come from laboratory tests of 
the inverse square law (ISL) such as the Eot-Wash experiment [14]. Since the bounds on chameleon theories coming 
from the Eot-Wash experiment have been discussed in great detail elsewhere [2, 3, 14] we have not focussed on them 
here, although it should be noted that many of the formulae derived in this work are also useful for estimating the 
size of any chameleonic signal that could be detected by that test. It should be noted, however, that the electrostatic 
shielding used by such tests, however, also acts as a chameleon force shield if the matter is coupled to the chameleon 
much more strongly than it is to gravity. Casimir force measurements and ISL tests therefore probe different values 
of the chameleon to matter coupling and taken together place the strongest constrains on chameleon fields with a 
gravitational or super-gravitational strength matter coupling. 

For m^^ ^ d < ti^^ and = or = V2 and A w Ao wo found the chameleon force between two parallel plates 
was generally about 0(1)% of the size of the total Casimir force at separations of d ~ 10 /zm. At larger separations, 
the dominant contribution to the Casimir force is expected to come from thermal effects the precise form of which 
is a matter of some controversy. Even still, in the parallel plate geometry, the chameleonic and total Casimir forces 
were found to be generally equal in magnitude for Add ~ 0(1) where A^ = Aq/A « 82 /um for A fa Aq. In the sphere 
plate geometry the situation is slightly more complicated but generally the chameleonic force dominates for smaller 
values of d. 

For A w Ao, the predicted chameleonic force per unit surface area is very small. It is therefore desirable to use 
larger rather than smaller test masses so as to increase the surface area and amplify any chameleonic force. It is no 
surprise then that the some of the best current constraints on it are provided by the 1997 Lamoreaux experiment [10], 
in which test masses with relatively large dimensions were used. The relativity large dimensions of the test masses 
also means that they have thin-shells for a larger range of couplings and so this experiment was sensitive to a larger 
range of coupling than were the tests reported in Refs. [13] and [17] where smaller test masses were used. 

We identified two future experiments with experiment with excellent prospects for detecting or ruling out chameleon 
fields associated with the energy scale of dark energy. Both experiments share some common features: they both make 
use of relatively large, 1 — 10cm, test masses, and both probe larger separations than previous tests. Additionally 
both experiments make us of torsion balances to measure the forces. The experiments in question the new tests 
proposed by Lambrecht et al. [26] and under construction in Grenoble and the new Lamoreaux experiment [25, 27]. 
We found that these tests are complementary to each other. Provided the laboratory vacuum is of sufficient quality. 
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the Grenoble experiment [26] has the sensitivity to detect or rule out almost all of the strongly coupled chameleon 
theories which are at the same time dark energy candidates (A = Aq = (2.4 ± 0.1) x lO^'^eV. However, before this 
could be done, the controversy that surrounds the correct method for calculating the thermal Casimir force would 
have to be settled, and the total Casimir force modelled theoretically to an accuracy of better than 1% at distances 
of w 5 — 10 fim. The Grenoble experiment has, in fact, be designed to help settle this controversy and so it is highly 
feasible that in the longer term they will be able to detect or rule out chameleon fields. 

The new experiment proposed by Lamoreaux [25, 27] could be used to probe even larger separations than the 
Grenoble test. For n ~ 0(1), A w Aq and V = Vi ot V = V2, we predicted that the chameleonic force would 
dominate over the total Casimir force when d > 33 /im. By probing separations slightly larger than this, the new test 
is predicted to have the sensitivity to detect chameleon fields with < n < 1.6 for a larger range of M than could the 
Grenoble experiment; although the Grenoble experiment sees a much larger range of M. Crucially though, since the 
chameleonic force would dominate over the Casimir force at these separations, the new test proposed by Lamoreaux 
could make such a detection without a detailed knowledge of the thermal contributions to the Casimir force. 

For m^^ ^ d < m'f^^ the chameleonic force that is detectable by Casimir force experiments depends is virtually 
independent of the strength with which the chameleon couples to matter, M, as well as the density and composition 
of the test masses. This, of course, means that whilst Casimir force experiments are excellent probes of the properties 
of the chameleon potential e.g. n, A and Aq, they are unable to provide much information about M. By a happy 
coincidence, however, the converse is true of experiments that search for the conversion of photons into chameleons, 
e.g. BRFT, PVLAS, BMV [9]. In the region of peak sensitivity the effects detectable by those tests depend only 
indirectly and as a result fairly weakly on V{(j)) but have a very strong dependence on M [6, 7]. Casimir force 
measurements and light propagation experiments therefore provide complementary constraints on the properties of 
chameleons. 

In summary: Tests of gravity generally assume that any new forces behave a lot like gravity: exhibiting a similar 
dependence on the density of test masses and being difficult to shield. The force mediated by chameleon fields obeys 
neither of these assumptions, and it is partly for this reason that chameleon fields arc so difBcult to detect or rule 
out using traditional tests of gravity. We have found that Casimir force measurements already provide the strongest 
constraints on chameleons fields with strong matter couplings. However, most chameleons theories where the potential 
is associated with the energy scale of dark energy, 2.4 x 10"^ eV, remain hidden from these tests. This underlines 
just how little we know about potential new physics associated with this small energy scale. The next generation of 
Casimir force tests however will offer greatly increased sensitive and will probe larger separations where the effect of 
chameleonic force relative to the Casimir force is more pronounced. Provided the controversy surromiding the thermal 
Casimir force can be settle, these experiments will have the sensitivity to detect or rule out almost all strongly coupled 
chameleon fields which are dark energy candidates. The prospect of detecting new physics, such as chameleon fields 
associated with the mcV scale of dark energy, in the laboratory is an exciting one. If no such detection is made then 
the potentials of viable chameleon theories would have to feature an energy scale <^ Aq which would seem highly 
unnatural. Whatever the outcome, by probing separations between 10 /im and 100 /im the new generation of Casimir 
force measurements have the potential to increase what is currently known about the nature and origins of dark 
energy. 
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APPENDIX A: GENERALIZED THIN-SHELL CONDITIONS 

In this appendix, we derive the condition for a body to have a thin-shell in a chameleon theory with general V{(j)). 
Consider an isolated spherical body with density radius i? in a background with density pi,. We take r be the 
distance from the centre of the body. The chameleon field evolves in the effective potential: 
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Wc have taken the body to be isolated. By this we mean that the length scale, Lf,, of the region where p = pb is large 
enough for 0(, far from the body. <j)b is defined to be the minimum of V^ff when p = pi,: 

V'{<l>b) = -^. (Al) 

We also define rub = m^{(j)b) = to be the mass of the chameleon near (/> = (pb- In the region where p = pb, 

perturbations in the chameleon field about (pb decay exponentially over a length scale ^ m^^. The isolated condition 
therefore translates to requiring uibLb ^ 1. 

The statement that a body of density pc has a thin-shell is equivalent to requiring that (p (pc deep inside the 
body, where V^g{^c', Pc) = defines pc i.e.: 

V'{4>c) = -^. (A2) 

We define rric = m^{(pc)- 

In this set-up, the chameleon field equation, Eq. (3), reduces to: 



d?(p 2d(p_ 
dr^ r dr 



2 r>J!, (A3) 



S + ?* = n«-n« r<R, (A4) 



dr^ r dr 
We define (po in r > Rhy the equation: 

\i (p K (pQ in r > R then the (p field produced by body behaves as if it were a linear perturbation about (p = (pb in 
r> R. 

We now prove that if ^ ~ (^o as r — > oo and < (pb then (p < (po in r > R. We take ^ = (pQ + ^i, and 
\4>il{4'o ~ 4'b)\ — > as r — > oo. We wish to prove that (pi < 0. We write (p\ = ^e~™'''"/r and then: 

[e-m,r^y ^ r{{V'{(P(, + (P,) - V'{(Pb) - ml{(Po + (Pi - (Pb)} • (A6) 

As r — > DO we certainly have (p < (pb- We have required V"'{(p) < which implies that ii (p < (pb then: 

V'{(P)-V'{(Pb)-ml{(P-(Pb)<0. 

It follows that as r ^ oo we must have: 

and so since e~"*''''x' ^ by as r ^ oo by the requirement that (/> ~ (/)o as r — > oo, we must have for large r that: 

■X^(r) > ^ x{r) < => (pi(r) < 

Thus (pi{r) < provided that: 

V'{(Po + <^i(r)) - V\(Pb) - ml{(Po + (Pi{r) -(Pb)<0 

but since required (pQ < (pb this certainly holds as r ^ oo and continues to do so provided that(jii(r) < 0. It follows 
(pi (r) < and thus (p < (po for all r > R. We will need this shortly. 

A full solution of the field equations requires (p{r ^ R) > (pc, since V'{(p) < 0. For a body to have a thin-shell, it 
must fulfil (p'^ (pc inside the body which requires that all perturbations in (p about (pc decay over a length scale that 
is smaller than R. A necessary condition for this is mcR > 1- It is enough, however, to have m{(p{r = R))R ^ 1 on 
the surface of the body. Since m{(p{R)) > m{(po{R)) (as (p{R) < (po by continuity) then its enough to predict that 
'm{(pQ{R))R ^ 1. If we assume that there is no thin-shell, then we would have to have m{(pf){R))R < 1, but if the 
assumption of no thin-shell leads us to predicting (po{R) < (pc then because of mcR ^ 1 this cannot be the case. 

We now find a suflacient condition for the existence of a thin-shell, by assuming that the body does not have a 
thin-shell and seeing when this leads to a contradiction. If a body does not have a thin-shell then inr > R: 

rr>„-mb(r-R) 

(P^(P^(r)^(Pb-— , (A7) 
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for some constant C. And in r < i? Eq. (A4) can be linearised to, 

dV 2 d(j) 



dr^ r dr 



mb{(t>-(t>b) 



Pb 



M 



and so in r < i? we would have: 



(j)b 



Mml \Mml J 



R sinli(m(,r) 
r sinli(m6-R) ' 



where 



e = 



Spce-'^"^ cosh{mbR) 



Mml 



= (F'(</.6) - V'(</),)) . 



-rriijR 



^ tanh(m(,ii) 

ITlbR 

cosh (mfcii) 



(A8) 



(A9) 



1 - 



tanh(mf,i?) 

ITibR 



and 5pc = Pc — Pb- However, if (^^{r = R) = 4>b — C 1^ 4'c then we found that there certainly is a thin-shell, which 
leads to a contradiction. A sufficient condition for the existence of a thin-shell is therefore: 



C>6b- 



(AlO) 



1. Linear Thin-Shells 



Consider Eq. (A8) which gives the form of </> in r < i? in those cases where it is acceptable to linearize the field 
equation about (p = (pb- If rribR ^ 1 then all variation of </> inside dies off in a thin region near the surface of the 
body (over thickness ~ 1/mb). Furthermore if this linearisation is correct then (j)c 4'b — Spc/Mm"^, and so ^ « 
deep inside the body. It follows from the definition of a thin-shell that such a body would have one. However, this 
type of thin-shell behaviour also exists in Yukawa scalar field theories where — const; it has nothing to do with 
the non-linear nature of the field equation which give the theory its chameleonic properties. We therefore deem this 
to be a linear thin-shell. 

If m,bR 3> 1 then a body will have a linear-thin shell provided (pc ~ (j>c which implies rUc ~ mfe. Outside of a body 
with a linear thin-shell, (p has the form; 

C^g-mfc(r-fl) 

<px<pb , 

r 

where C is given by Eq. (A9). Since C depends on 5pc/M it is clear that the chameleonic field far from the body 
depends on both the density of the body and the strength with which the chameleon couples to it. This is precisely 
what one naively expects to see in scalar field theories with a coupling to matter. This behaviour is, however, very 
much associated with linear field equations and it is not, as we shall see, what one finds when the non-linear nature 
of the chameleon field equation dominates the behaviour of the field. 



2. Non-Linear Thin-Shells 



In addition to linear thin-shells, which exist even in theories with a simple Yukawa scalar field theories, chameleon 
theories also exhibit non-linear thin-shell behaviour. If a body has a thin-shell of any description then all variation 
in (p is exponentially small outside a thin region near the surface of the body. As was shown in ref. [3] for V oc 
theories, when a body has a non-linear thin-shell, the chameleon field far from the body is independent of 5pc/M. 
This is one of the key features of chameleon theories that allows models with M <C Mpi to evade the constraints 
coming from experimental tests of gravity and searches for WEP violation [3]. 

Outside the body, (p obeys Eq. (A3) and by integrating this once we find: 

K^)' " V i<f')-v (<!»>) -V'{<l>b){<l>-<l>b) + J^^^(^^Jdx, (All) 
In r > i?, we have d(f)/dr > and so defining d = r — R, 
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and Us = y{r = R) we have: 



where 



Defining 



it is useful to write: 



X = I W{z)dz, 



W{y) = {V^-V' + 2V2y/r)-\ 



(A12) 



-f 



2 fd4> 
dx 



dx. 



W{y)-^= -V- mb^/2{y^ - x^) 



+ 



rribV'^iy'' - X') + 



2V2y' 



(A13) 



If V = Vb + V^{(j) — (f>b) + \m1{(j) — (f>bY then the terms in the first set of square brackets would vanish (note that 
<j)b> (j))] these terms are therefore associated with the non-linear nature of the field equations. Whatever the potential 
is, the terms in the second set of square brackets arc manifestly increasing more slowly than y as y — > oo. From the 
definition of y, it is clear that if near the surface of the body: 



vi - y » 



2^/2y _ 2 d(^ 
r r dr ' 



then <j) is varying over scales that are small compared with r, and so we have thin-shell behaviour of some description. 
If additionally: 



VI -V- m6 V2(t/2 - X') 



m6V2(y' - X') + 



2V2y 



(A14) 



This is 



near the surface of the body then the non-linear terms in the field equations dominate the behaviour of 
therefore the condition for a non-linear thin-shell. 

If the non-linear thin-shell condition holds then <f) is a quickly varying function near r = R. Let us define 5R ^ R 
to be the scale over which d(/)/dr varies near r = R. We then have ~ 0{Ar/R)y^ 
of a body for which Eq. (A14) holds. For r <^ Rwe then have: 



and so 2/^ near the surface 



-V- mbsJW 

2 -2 

y «y 



X") « Q{y) = vi -V - mbV2j/, 

= v{4>)-v{ct>b)-v'{<i>b){ci,-4>b), 



(A15) 
(A16) 



which equality as y 



oo. It can then be checked that V" < implies that: 

dlnQ 



dlny 



>1, 



and so Q increases always increases faster than y. Thus when Eq. (A14) holds, W{y) — > faster than 1/y as y 
increases. Thus not only docs the integral in Eq. (A12) converges as yg oo, but this should additionally provide 
a very good approximation to y{x) whenever y{x) <C ys and Eq. (A14) holds for y = y^. The far chameleon field 
perturbation produces by bodies for which Eq. (A14) holds near their surface is therefore almost independent of t/g. It 
is for this reason that the perturbation in the chameleon field far from a body with a non-linear thin-shell is virtually 
independent of the value of (f> on the surface of the body and, as a result, of Pc/M. 

We now return to the non-linear thin shell condition, Eq. (A14), and rephrase it in a more useful form. We begin 
by assuming that a body does not have a thin-shell and hence that (j) ~ </>o(»') inr > R. We then look to see when this 
assumption leads to a contradiction. The assumption that there is no thin-shell and cj) ^ 4>o implies, by Eq. (A14), 
that at r = i?: 



V'{<Pb)-V'{(j>b-C) 



< 



2C_ 
R^ 



{l+TUbR- 



mlR^) 



(A17) 
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Using Eq. (A9) to give C for a non-thin shelled body, we find that Eq. (A17) is violated, implying that there must 
be a non-linear thin-shell, when: 

where 



/(mfci?) = 2e-™''^cosh(m6i?) 



mbR mlE? 



tanh 



mbR\ 



rubR ) 



If m\,R <C 1, f{mbR) ~ | and if rubR ^ 1 then /{nibR) ~ 1; finib) is monotonic in mbR. The non-linear thin-shell 
condition (or the thin-shell condition for short) is therefore approximately equivalent to: 

^ _ (Pc - Pb).f{mbR)R^ > (f)b-(p - ~ fi^bR)) ^^^g^ 



or 



2M[mlR^ + mbR -M] ~ ^ ^ Mm^ 



mlR^ ^ 2Mml {ct>b - ^c) _ 2(1 - f{mbR)) ^ ^ ^^^^^ 



mjR^ + mbR + 1 ~ (pc - Pb)f{mbR) f{mbR) 

These conditions could, alternatively, be written as C > Cthin where 

Ctbin{l + mbR + ml R^) ^ ^ 
V'{cl>b)-V'{cj>b-Cti,in) 2- 

It is clear from its definition that the Cthin depends only on R, mb and the form of V{(j)). It is independent of Pc/M. 

Since, in a thin-shelled body, almost all variation in (f> takes places in a thin region near the surface of the body, the 
non-linear terms in the field equation should only be important over scales <C R. Far from a thin-shelled body, 
should therefore behave if the field equations were linear i.e.: 

^ ~ ^6 , 

r 

for some C*. If C* <C Cthin, then the thin-shcU conditions would not be satisfied near the surface of the body, and 
if C* Cthin then by Eq. (A14) the non-linear terms in the field equation would begin to dominate far at r -R. 
It must therefore be the case that C* « Cthin for thin-shelled bodies. As r/ii — > oo then we have, for a thin-shelled 
body, that: 

cP^cPb-^^^ . (A20) 

r 

As should be the case, the large r behaviour of R is, to leading order, independent of pc/M and depends only on mb, 
R and, through Cthin on the form of V{(l)). 

APPENDIX B: FORCE BETWEEN TWO PARALLEL PLATES 

It was shown in section IV A that the chameleonic force between two parallel plates with thin-shells and separation 
d is given by: 

^ = viM - v{ci>b) + v'{<Pb){<l>b - M < viM - v{^b), (Bi) 

where (provided d 1/mc, where mc is the chameleon mass inside the plates) (j)o{d) is given by: 

f°° d 
^/2j^ W{y;ct>o)dy=-, (B2) 



where we have defined y = ^/V - Vq - Vl{(l) - M and l/W{y; 0o) = {V^ - 1^' (</>)) > 0. We define P{y) by: 

d\nW ^ 2V"{V-V{M-Vl{<j>-^o)) ^ _ 
d\ny (y/(^)_v,')2 - ^y'' 
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As y ^ oo, the condition that V" < 0, implies that P{y) > 1 as y ^ oo. Wc define y^i by P(y_i) = 1; the dominant 
contribution to the integral in Eq. (B2) then comes from 0{y-i) values of y. Defining by y{4>-i) = y-i, we find 
a good estimate of <p-i by expanding P{y) about y to order — ^o)^: 

Fo'>_i-^o)'«yo'-V^6- (B3) 
Now, linearization of the field equations about (po is a good approximation provided that: 



«: 1. 



The point (p — (p-i therefore lies close to the edge of the region where the linear approximation is valid. Since this is 
the case, we approximate the integral in Eq. (B2) by expanding W{y) out about y = = ^q): 

- V'f = 1/W\y) ~ (n' - V^f{l + 2aV + k'^aY + 0{y^)). (B4) 

where we have defined 



If A;^ < 2 then to 0(y*) we can rewrite 1/W'^{y) as: 

« - [(1 + + Oiy')] , (B5) 

where = a^/p and {2p -l)/p = fc^, and < fc^ < 2 implies p > 1/2. 

Wc define Ueff = (2 — k'^)/{k^ — 1) so that p = 1 + l/ncff. W{y) drops off faster than 1/y for y^ > l/k"^. 

The approximation of W{y) given by Eq. (B5) is therefore approximately the same as that given by Eq. (B4) for 
y'^ < as long as the new 0{y^) terms introduced in Eq.(B5) are smaller than the 0{y^) terms at y^ = 

This requires: p > 4/7 i.e. k'^ > 1/4. To ensure that we are likely to be justified in ignoring the 0{y^) terms, we 
therefore require k^ > 1/3. Therefore for 1/3 < fc^ < 2 we find that: 

d V2 f°° dy 



2 V^-V^Jo (l + c2y2)i+i/ne«- 

Performing this integral, we arrive at: 



rieff \2 2 Ueff 

where B{-, •) is the Beta function, Ucs = (2 - A:^)/(fc^ - 1) and 1/4 < fc^ < 2 i.e. ricrr > or neff < -7/3. When 
y = + A^(A/(/))" and toq ^ nib this approximation for mod is actually exact and rieff = n; note that for such a 
potential we require n < — 4 or n > for a valid chameleon theory to emerge. Note that the condition mo » rrib 
implies: 



y neff V2 2 rieff/ 

We now consider the small fc^ case. If fe^ is small it is either because fcg = V^'V^/niQ is small, or because 1 — V^Vq 
is small. In the former case, the approximation: 

^ (yj - Vlfil + 2aV + k'aV + 0{y')), 



Wiy) 



is generally no longer valid to 0(t/*) at the point y = where W{y) ^ C/y, for some C. When this approximation 
does hold at y = y-i, we have a?y^i ~ l/fc^. The O(y^) terms are of the order —gk'^a^'y^ /3 where, 

.(4) 
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For the 0{]f) term to be smaller than the 0{]/) at y = ?y_i we need l.gfc^a^y^ J = \g\ < 3^2 y{'i)y^,^y(3)y-(2) ^ 
0(1), then this condition generally holds provided that > 1. It is for this reason that we took k-^ > 1/3 in the 
previous case. If fc^ is smaller than a 1/3, however, then it is generally the case that \g\ < 3k'^ does not hold. If fc^ 
is small because k^ is, then the behaviour of W{y) near the point where it drops off like l/y is determined by higher 

derivatives in V than yj'*'' . However, chameleon potentials with k^ <^ 1 are not generally very natural, and so we do 
not consider them here. What is important to understand is behaviour of (f)o as ^ ^b- 

As d ^ 00, (f)Q converges to cf)b and so fc^ 0. It is important to understand how (po and hence F^jA behave when 
d > 0{l/mb). For large separations, 0o 06- For close to 4>o then we have: 

y' + yl^%cp-<i>bf + ^{<t>-<i>bf, 



where 



Expanding W^{y) we find: 

1 



We therefore have: 



yo = Vo-Vb-V^{cPo-cPb) = ^. 



= {Vb - Vy « mt{ct> - 4>bY + mtvricp - <t>bY + 0{{<t> - ^bT). 



2ml{y' + yl)-'-^{2{y^ + yl)fl\ (B7) 



W^{y) 3mb 
We then have: 

nibd f°° dt 



2 Jo V^Ti^/i + iWTiy^' 

where 

3ml ~ 3 

For rubd 1 we then have: 

rubd f°° dt 



f 

Jo 



VT+7^"cosht 



2sinh-i ^ - ^ + 0{l^) « ln(8) - 21n(Z) + 0(/*). 



Thus as d 00 we have: 

nibd 



ln(12) - ln(A;^). (B8) 



This approximation therefore certainly requires < 1 to be valid i.e. fc^ < 3/2. Furthermore we must also have 
mb « toq. If this holds then: ttiq « ml{l + k^) and so mg/mj « 1 requires: fc^ <C 1, which in turn requires 
nibd ^ 21n(12) « 5. Thus for mjrf 5 we have 



Finally, we consider the behaviour of mod and hence F^/A for k'^ > 2. In this case we write: 

1 



(B9) 



W^{y) 



(Vo' - + 2a'y^ + 2aY + {k' - '^cYy' + 0(y')), 
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and so to 0{y'^) we have: 



We therefore find that: 



mod 



W^{y) 



r3/2 



{Vi-Vife^'^'y\l + {e-2)aY). 



/2 



/4 



2^2(^2 -2)( 1/2) 

where J_i/4 and y_i/4 are Bessel functions. For small 4(A;^ — 2) this gives: 

mod K \/2^ (1 - 3(A:^ - 2)/8) , 

and if fc2 ^ 2 we have: 

5.24 



/4 



mod i 



^4' A) 



(BIO) 



where, as above, -B(-, •) is the Beta function. Note that the 1/3 < /c^ < 2 and k"^ > 2 approximations for mod{k^) are 
continuous at A;^ = 2. 



APPENDIX C: FORCE BETWEEN A SPHERE AND A PLATE FOR POWER-LAW POTENTIALS 

WITH < n < 2 

In this appendix we derive the force between a sphere and a plate for a power-law potential with < n < 2. The 
potential is taken to have the form: V = + A''(A/^)". When mo '> nib, Eq. (22) is exact and we have: 

l2n + 1 /I 1 l\ n + 2 

mod = \ B -,- + -= ^n{n + l)i^„=" . 

V n \1 1 n J 

This result is valid so long as: 

/2n + 2„ /I 1 1\ 

For < n < 2 the RHS of the above expression takes values between about 2.5 and 3.2. Provided mbd <C 3, the force 
dF^ between two parallel surfaces with separation s and area AA is: 

^«if„A^(A5)-^. (CI) 

Let us now take one of the surfaces to be a sphere of radius R. The minimum separation of the sphere and plate 
is taken to to be d. Provided mbR <C 1, <C -R, we showed in section IV B that the total force, -F|°*, is, to a good 
approximation, given by: 

Ff'^2^Rj^^^ds. (C2) 

When n > 2 and toq S> mf,, dLF^{s) / dA drops off faster than 1/s. No matter what n is, wc found in appendix B 
that when toq ~ m;,, dF^[s)/dA oc e^™''''. Therefore when n > 2, the dominant contribution to the comes from 
values of s « d. If, as is often the case, rribd -C 1, we can therefore calculate to leading order knowing only the 
form of dFff,{s)/dA for rribS <C 1. This is given by Eq. (CI). This calculation is performed in section IVB and so is 
not repeated here. 

In this appendix we consider the behaviour of F^ when < n < 2. In these theories, when rribS <C 1, dF^{s)/dA 

drops off more slowly than 1/s. When rat,s > 5, dF(f,{s)/dA oc e~™'''^ . Wc therefore expect that F^{s)/A first drops 
off faster than 1/s when s = s_i ~ (^(l/mb) < b/mb- The dominant contribution to the integral in Eq. (C2) then 
comes from values of s that are 0{s-\) i.e. 0{l/mb). 
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We found in Appendix B that the assumption, toq ~ mt, required rribd > 5. It is therefore safe to assume that 
mo(s_i) is large compared to m^. The 1/3 < < 2 approximation for rngs, Eq. (B6) can therefore be apphed. We 
define x = V^/Vq; a; — > 1 as s — > oo and a; — > as s ^ 0. We then have: 

n+ [n + ■2):i: 

neS = -. -, — — ;r;— ) 

1 — (n + 2)x 

and so by Eq. (B6) we have: 



m,s = x-^mos « ^ x^ J ^^P^B f 1, -1 + " + ) . (C3) 

]^n+{n + 2)x \2' 2 n+{n + 2)xj ^ ' 

Now 

^ = Vb - H - yi{<^^ - <^h) = H (ar-* - 1 + n(x^ - 1)) . (C4) 

By Eq. (C2) 

Fr=2.R r'jm,s-2.R r^d. 

* Jo dA Jo dA 

For < n < 2, and provided nibd <C 1, the second term on the right hand side of the above expression evaluates to: 



Jo dA 2 — n 



We define: 



Using Eq. (C4) we have: 



Jo 



Fo = 2nRl ^ds. 

dA 



Fo = 2ttR— [ - 1 +n(a;^ - 1)) d(mbs), 

mb Jo \ ' 

= 2itR— [ (x-^ - 1 + nix^ - 1)) L'(x)dx, (C5) 
rrib Jo ^ ' 

where the last equality follows from Eq. (C3) and L'(x) = dL/dx; L{x) is given by Eq. (C3). 

We now wish to evaluate Fq approximately. The dominant contribution to Fq comes from values of x near the point 
where: 

(x-^ -1 + n(a;^rTT _ i)^ L{x) (C6) 

takes its maximum value. For all < n < 2, this maximum occurs for values of x < 0.5. Our assumption that 

mb/ms = a;("+2)/(2("+i)) <|; 1 is therefore justified. 

We define aimax to be the value of x at the maximum of Eq. (C6); the smaller n is, the larger Xmax becomes. For 
< a; < 0(a;max) we find that: 



n+(n + 2)a; V2 2 n+(n + 2)a;y V n V2 2 nj l-x 



where 



3 =^ 
^" 2n2 



2(n+l)(*(i)-*(l + i 



and where 5" is the Digamma function i.e. ^{x) = T'(x)/V{x). We use Eq. (C7) to approximate the right hand side 
of Eq. (C3), and thus to approximate L(x) as: 



2' 2 nJ (l-x) 



38 



Putting this approximation for L{x) into Eq. (C5) we find 



where: 



_ Mn+1) / 
^" (n + 4)(n + 2)V 

We note that D2 = 1. Fq may therefore be written as: 



2- n 
3(n + 2) 



Fo«27r 



n + 2 
2-n 



rub 



where we have defined 



n V2 2 n 



Thus for < n < 2 the total force between a sphere and a plate with separation d <C m^, ^ <C i? is: 

I- 9 — «. -1 



J^^°*(0 < n< 2) = 27rA^ RKn 
By taking the limit of this equation as n ^ 2 we find if n = 2: 



- (Arf)^ 



Ff\n = 2) = 2nK^RK2 



log 



where D'^ = d£)„/dn|„=2. We find that AD'^ = 2 + 21og2, 02 = 2\/3 and K2 = \/2. Thus: 

Fr(n = 2) = 2V2.A3i?log(^). 



(C8) 



(C9) 
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